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THE SYSTEM K.0-Mg0-SiO;' 
PART 1 
EDWIN WOODS ROEDDER 


ABSTRACT. The phase diagram of the more gc:Jogically significant 
portions of the condensed system has been determined by the quenching 
technique. Four new ternary compounds were foun: and their thermal re- 
lationships determined. Their compositions are K,0-5Mg0O.12Si0,, K,O. 
MgO.5Si0O,, K,O0.MgO-3SiO,, and K,0-MgO-SiO,(?). The first of these 
may occur naturally but has not been reported, possibly because of its 
optical similarity to quartz, The second and third appear to be isomorphous 
with leucite and kalsilite, respectively, and the last is of comparatively 
little geologic interest. Thermal, optical, and X-ray data are given for 
these compounds, and composition and temperature have been determined 
for 28 of the invariant points occurring in the system, Experimental 


technique, crystallization behavior, and geological application are discussed 
briefly. 


INTRODUCTION 

HE immediate purpose of this study was to aid the investi- 

gation of the four-component system K,O-Mg0O-Al,0,-SiO, 
by Dr. J. F. Schairer,’ which was studied because of its bearing 
on the stability relationships of the magnesian minerals for- 
sterite, clinoenstatite, and cordierite, and the potash-alumina- 
silica minerals leucite, potash feldspar, and quartz. Four 
ternary systems limit the quaternary system K,O-MgO-Al,0,- 
SiO,. One of these, K,0-MgO-Al,0;, has not been investigated 
because of the lack of direct geologic interest; two others, 
Mg0O-Al,0,-SiO, (Rankin and Merwin, 1918) and K,0-Al,0,- 
SiO, (Schairer and Bowen, 1947a), have been studied previ- 
ously ; the fourth, K,0-MgO-Si0., is the subject of this report. 
Although there is no naturally occurring mineral known con- 
sisting of K,0, MgO, and SiO,, during the course of the work 
several new ternary compounds of interest to petrology and 
mineralogy were found and are reported here. 


1A dissertation submitted in partial fulfillment of the requirements for 


the degree of Doctor of Philosophy in the Faculty of Pure Science, 
Columbia University, 1950. 


2 Preliminary report in abstracts of papers presented at the meeting of 
the Geological Society of America in New York, 1948. 
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Potash, magnesia, and silica are among the seven most 
abundant constituents of igneous rocks, and Washington 
(1915) has shown a direct relationship between the amounts 
of potash and magnesia in certain rocks. The system is of 
interest also in the broader problem of the chemistry of silicate 
formation and stability. Only that portion of the condensed 
system higher in silica than the join K,0°Si0,-MgO was 
studied. 

In addition, the work was part of a long term project being 
carried on at the Geophysical Laboratory on the stability re- 
lations of the rock-forming dark micas—the biotites (Adams, 
1947). As biotites are essentially hydrous, potash-magnesia- 
alumina silicates with a partial substitution of Fe** for 
Mg* *, their compositions lie within the six-component system 
K,0-Mg0-FeO-Al,0,-SiO,-H,0. Since the presence of water 
and iron oxide causes considerable experimental difficulty, the 
most reasonable approach was a study of the system K,O- 
MgO-Al,0,-SiO,. A knowledge of the limiting system K,O- 
Mg0O-Si0, is essential to an understanding of this quaternary 
system. The only previous work on the system K,O-MgO-SiO, 
involved the addition of water as the fourth component to a 
small range of compositions with the formation of hydrated 
mica-like clay minerals (Noll, 1944). Another possible ap- 
plication of the system K,0-MgO-SiO, is in the further study 
of the system MgO-SiO,-H,O (Bowen and Tuttle, 1949) 
where the addition of K,O may aid the formation of a liquid 
phase (Adams, 1947). 

This work was done at the Geophysical Laboratory of the 
Carnegie Institution of Washington on a one-year fellowship. 
The writer wishes to thank the Director of the Laboratory, 
Dr. L. H. Adams, for this opportunity, and Dr. J. F. Schairer, 
under whose direction the work was done, for his aid and 
encouragement. The writer is indebted also to the other mem- 
bers of the Laboratory for many helpful suggestions, and 
particularly to Drs. Schairer, Morey, and Yoder for a critical 
reading of the manuscript; however, the writer alone is re- 
sponsible for any errors which may occur. He also expresses 
his appreciation for the cooperation, assistance, and advice 
of the members of the staff of the Geology Department of 
Columbia University, particularly to Professor C. H. Behre, 
Jr., who initiated the original discussions leading to the au- 
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thor’s connection with the Geophysical Laboratory; to Pro- 
fessor Paul F. Kerr for his cooperation in making this arrange- 
ment possible; and to Professor Ralph J. Holmes for advice on 
X-ray diffraction techniques. 


EXPERIMENTAL PROCEDURE 


Raw Materials——The starting materials for most of the 
compositions were KHCO,, quartz, and MgO. The KHCO, was 
of reagent grade,® and was prepared by crushing to a powder 
and drying over KOH. The quartz used was the regular Geo- 
physical Laboratory supply of Lake Toxaway quartz, crushed 
and leached with HCl. Analysis by J. W. Greig gave a residue 
of 0.04 per cent (using 5 grams) on evaporation with HF 
and H,S0Q,. It was prepared for weighing by heating for two or 
more hours at 1500° C.* This heating serves to eliminate liquid 
inclusions and converts the quartz to cristobalite, facilitating 
subsequent grinding in an agate mortar. In a few special 
batches Mallinckrodt “Special Bulky” silicic acid was used 
because of its fine grain size,’ but in these cases a water analysis 
had to be made simultaneously because of its large and vari- 
able amount (15.5 per cent). The MgO used was reagent 
grade,® although analysis by J. F. Schairer showed 0.2 per 


8J. T, Baker Chemical Co., Lot no. 12738 analysis: 


Heavy metals (as Pb) 
Ca, Mg, and NH,OH ppt. 
+All temperatures are in degrees Centigrade. 
5 The powder has a bulk density of only 0.11 grams per cc. 
6J. T. Baker Chemical Co., Lot no. 82735 analysis: 
per cent 
Insol. in HCl 0.02 
Loss on ign 25 (1035 E. R.) 
Soluble in H,O 0.4 
Heavy metals (as Pb) 
NH,OH ppt. 
(NH,),S metals (as Fe) 
Sulfate and sulfite 
SiO, 
Cl 
NO, 0.005 
Ca 0.05 (<0.2 J. F. 8.) 
Ba 0.005 


per cent ra 

cl 0.008 

PO, 0.001 

NH, 0.000 

Fe 0.000 

SO, 0.001 
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cent total alkalies. It was prepared for weighing by heating 
to 1500° for one hour. This minimizes water adsorption, but 
corrections still had to be applied to each portion weighed, 
as determined by individual weight-loss measurements upon 
ignition of each portion weighed. 


Preparation of Potassiwm Silicates—In almost every case 
the KHCO, was first converted to a pure, completely crystal- 
line potassium silicate compound (K,0°4Si0,, K,0°2Si0,, or 
K.,0°SiO,) before using in order to minimize volatilization of 
K,O during the preparation of the various compositions. A 
glass of the composition of the desired potassium silicate 
compound is first made by mixing together on low-ash filter 
paper the stoichiometric amounts of KHCO, and SiO, to give 
20 grams of the silicate. The mixture and the paper are put 
into a 40 cc. platinum crucible, carefu.ry heated up to 1000° 
over a predetermined 48-hour heating schedule to minimize 
foaming and volatilization of K,0 (Kracek, 1932, p. 2529), 
and finally heating for two hours at 1300°. When cool it is 
weighed ; any loss in weight in excess of the CO, and H,0 given 
off by the KHCO, is attributable only to volatilization of K,0 
and must be corrected for by the addition of the calculated 
amount KHCO,.' Furthermore, the batch is far from homo- 
geneous at this stage and must be carefully removed from the 
crucible, crushed, replaced, and remelted again for two hours 
at 1300°, the process being repeated until a completely homo- 
geneous glass is obtainc? he extremely hygroscopic nature 
of this glass and the persistence of small amounts of CO, 
make is necessary to crystallize it, this being done by holding 
it at appropriate lower temperatures with occasional crush- 
ings until it remains a completely crystalline powder at the 
temperature of the furnace. As long as uncrystallized glass 
is present the powder will sinter together in the furnace. This 
crystallization is completed in 3 to 4 days,® after which the 


7 The amount of this volatilization of K,O is small—usually less than 
10-20 mg. in a 20 gm. batch—but sufficiently variable due to irregularities 
in the process that it must be corrected for in each batch separately. After 
this first fusion at 1800° further volatilization is found to be insignificant, 
and the volatility of silica at these temperatures is practically nil (see 
Sosman, 1927, p. 109). In those cases where the metasilicate, K,O-SiO,, 
was made up, failure to eliminate all of the CO, caused troubles which 
will be considered beyond under “Alkali Losses.” 

8 “Seeding” with previously obtained crystals has to be resorted to 
in the case of K,0-4SiO, because of its sluggish crystallization. 


4, 
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crystalline powder may be kept over KOH with little gain in 
weight except in the case of K,0'SiO,. 


Preparation of Batches.—All except the high alkali com- 
positions were determined entirely by synthesis. The calculated 
amounts of potassium silicate, magnesia,’ and silica were 
weighed out to make the individual compositions listed in table 
1. Ten grams was a convenient amount to make and should be 
less affected by systematic errors. Each batch was carefully 
mixed on low-ash filter paper and heated to 1300-1500° in a 
15 cc. platinum crucible for one or more hours,’® depending 
upon the viscosity of the liquid, then quenched, crushed,” re- 
placed in the same crucible, and fused again. This process was 
repeated until a homogeneous glass was obtained. The homo- 
geneity was checked by microscopical examination of the 
crushed glass when immersed in an index oil exactly matching 
the glass, making use of the more rapid change of the index of 
refraction of the oil with change in temperature to obtain a 
perfect match. Extremely small local variations in the com- 
position of the glass result in index variations which are 
rather prominent by this method. Viscous batches may take 
five or even more such fusions alternating with crushings be- 
fore becoming homogeneous.”* As has been emphasized else- 

® The small amount (0.05-0.10 per cent) of water taken up by the MgO 
during storage over KOH (following the 1500° ignition) was corrected for 
by weighing out the magnesia first into the crucible, igniting at 1300° for 


15 minutes, reweighing and correcting for the loss in weight before adding 
the other components. 


10 Weight loss determinations on a number of such batches (using a 
smooth plate cut from a corundum crystal as a plaque) showed practically 
neglible volatilization of alkalies, as changes in weight were within the 
limits set by the vaporization and condensation of platinum in the furnace 
at such temperatures. In the high alkali portion of the system weight 
losses were sizable (see “Alkali Losses” beyond). 

11 The laborious process of fine grinding of many samples is made 
considerably easier through the use of a Plattner-type hardened steel 
“diamond” mortar with the crushing end of the pestle of smaller diameter 
than the main pestle shaft and eccentric to it. The use of a steel mortar 
is recommended in place of other materials because contamination of the 
sample can be kept to a minimum through magnetic separation after each 
grinding. It has been shown (Sandell, 1947) that, without magnetic 
separation, grinding of quarts in a Plattner-type mortar yields 0.02 to 
0.03 per cent Fe in the ground material. 

12 The higher alkali glasses cannot be made perfectly homogeneous by 
this method as small variations in index persist indefinitely, presumably 
due to volatilization of alkali from the surface. The later crystallization 
of these glasses, however, forming crystals with individual compositions 


different from that of the glass, represents a far greater degree of 
“inhomogeneity” than occurs in the glass. 
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where, this homogenization is an important and essential step 
if true equilibrium is to be obtained in the subsequent quenching 
experiments with the material. 


Crystallization of the Batches.—All compositions were pre- 
crystallized partially or completely before quench runs were 
made.** To do this a portion of the homogeneous glass obtained 
above was held at temperatures which were 50 to 100°** below 
the approximated liquidus (as predicted by experience and 
extrapolation from nearby points). The viscosity of the more 
viscous glasses increases rapidly with decrease in temperature, 
hence the crystallization velocity is so small as to be experi- 
mentally impractical if crystallization is attempted at much 
more than 100° below the actual liquidus, and, of course, no 
crystallization can occur if held above the liquidus tempera- 
ture. Even under optimum conditions (correct temperature, 
frequent crushings, etc.) batches in the vicinity of the two 
lowest ternary eutectics did not begin to show visible crystals 
(1-2 microns in length) within the first month of crystalliza- 
tion. Such batches could not be completely crystallized in the 
time available for these studies, but many of the others could 
be, and these completely crystalline preparations were useful 
in determining eutectic temperatures. The viscosities of the 
glasses in a few cases were found to be so high that the 
powdered glass would not flow together to a coherent mass 
even when held for 24 hours at temperatures near the liquidus. 
Elsewhere in the system the viscosities were lower and liquidus 
temperatures were higher, so crystallization was usually initi- 
ated within a few days or even within a few seconds (in the 
very high temperature, low silica melts). Microscopical ex- 
amination at frequent intervals during the crystallization 
provides a very rough but useful indication of the phases and 
temperatures to be expected during the quenching experi- 
ments. Wherever possible coarse crystallization was avoided 

13 In most silicate studies this is essential because of the sluggishness 
of crystallization compared with melting. In some cases it was found that 
compositions which required months at optimum temperatures to form 
the first visible crystal nuclei and weeks more, with frequent crushings, 
to convert even a small percentage of the glass to crystals, would attain 
equilibrium in less than a week in the quenching runs, where equilibrium 
is approached by melting rather than by crystallization. 

14H. R. Swift (1947) found that 40-100° gave the highest rate of crystal 


growth for a number of silicate compounds but that there was wide varia- 
tion. See also the discussion of crystallization velocity in Tammann (1925). 


| 
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because of the long time necessary to attain equilibrium be- 
tween crystals and liquid during the quenching experiments 
with such material. A few of the more viscous melts with very 
low liquidus temperatures had to be subjected to “hydro- 
thermal” conditions*® in order to induce crystallization, but 
in these cases subsequent “dry” crystallization was usually 
possible (on another portion of the glass) after the cor- 
rect temperature range had been determined on the “hydro- 
thermally” formed crystals. 


Quenching Technique.—All of the equilibrium data was ob- 
tained by the quenching technique (Shepherd, Rankin, and 
Wright, 1909; Ferguson and Merwin, 1919; Faust, 1936; 
and Hill, Faust, and Reynolds, 1944). A charge of about 10 
mg. of the crystalline material in a small 0.01 mm. platinum- 
foil envelope’* is suspended in the hottest portion of the furnace 
(region of minimum thermal gradients), with a thermo- 
couple junction about 1 mm. to the side of the center of the 
envelope, for a length of time sufficient to attain equilibrium 


15 The method used involved sealing two cylindrical loosely-covered, plati- 
num-foil capsules into a fused silica tube. One capsule contained approxi- 
mately 0.2 gm. of the powdered glass to be crystallized and the other about 
the same amount of ground obsidian. The obsidian used was that from 
Mono Craters, California, containing approximately 0.07 per cent F, 0.13 
per cent Cl, 0.2 per cent H,O*, an unknown amount of H,S, and with SiO, 
and Al,O, 76 and 12 per cent respectively (see Shepherd, 1938). A sample 
of the ground obsidian was fused and showed a 0.41 per cent weight loss 
when ignited over a Meker burner. 

The whole tube was then put into the furnace at about 100° below the 
estimated liquidus of the glass to be crystallized. The heat drives some of 
the volatile “mineralizers” out of the obsidian and into the atmosphere of 
the tube, from which a portion is absorbed by the glass: sample. Although 
the partition of these volatiles between the obsidian, the tube atmosphere, 
and the glass sample is not known, the glass can hardly have more than 
one-third of the total amount present, and may possibly have much less, 
but this small amount of “volatiles” dissolved in the glass is adequate to 
induce crystallization in extremely viscous glasses in 24 hours, where 
holding for one month at the correct temperature in the “dry” condition 
(ie., with only the partial pressure of water naturally occurring in the 
air in the furnace) caused no crystallization. 


16 To make these envelopes a piece of foil approximately 10 x 16 mm. 
is folded so that the 10 mm. edges are together. The 8 mm. sides are 
then folded back over themselves twice, taking about 1 mm. each time. The 
resultant 6 x 8 mm. envelope is spread open at the top, the charge is 
inserted, and the top is then folded down. A pinhole through this top 
fold permits hanging. If the folds are pressed tightly, they will not leak, 


except in the case of high-alkali compositions where a double envelope 
is needed. 
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and then dropped out of the furnace into cold mercury by 
suddenly melting a thin (0.2 mm.), supporting, platinum wire 
by an electric current. The “quenched” foil envelope is torn 
apewt and the crushed sample examined with standard petro- 
graphic immersion methods (and X-ray where necessary and 
possible) to determine what crystalline phases were present 
at the temperature of the furnace, any liquid present having 
quenched to a glass. In a few cases, particularly at high 
temperatures and where quenching was not as rapid as de- 
sired, feathery crystals formed around each of the euhedral 
crystals present in the glass, or around other nuclei such as 
the thin platinum crystals that form on the surface of the 
melt from furnace vapors. Because the amount of such crystal- 
lization varied directly with the quenching time, the feathery 
crystals obviously formed on quenching and could be ignored. 

One minor change in the usual technique permitted the heat- 
ing of a larger number of different quench charges simultane- 
ously within the same furnace.’* Previously porcelain rings 
were used as separators between the platinum envelopes of such 
charges to prevent the welding together of the platinum sur- 
faces. The substitution of slips of ordinary white glazed paper 
for these porcelain rings was found satisfactory, the ash from 
the paper being sufficient to keep the platinum surfaces from 
touching at temperatures below 1040°. Above this tempera- 
ture fusion of the paper ash makes it necessary to substitute 
filter paper slips, dipped in a thin suspension of MgO. With 
this procedure 10 runs could be made simultaneously whenever 
desired. Even with 10 envelopes, the resultant bundle is smaller 
than the former maximum bundle of 4 envelopes and 3 porce- 
lain rings, and also more thermally conductive, resulting in 
smaller thermal gradients throughout the bundle and more 
rapid quenching. 


Determination of Equilibriwm.—The attainment of equilibri- 
um in the quench runs was established by several criteria: 

1. Absence of rounded or corroded crystals—this is not a 
reliable criterion either as positive or as negative evidence but 
is a useful indication. 

11 Whenever possible several different charges were run simultaneously 
to save furnace time. The bundle of charges was sufficiently small to be 
entirely within the region of neglible temperature variation, a cylindrical 
space approximately 1 x 1 cm. By this procedure differences between 


batches could be detected even though they were smaller than the accuracy 
of the calibration of the thermocouple. 


} 
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2. Agreement of data with nearby points and with other 
data on the same point—this is particularly useful. 

3. Duplication of results by starting with the same com- 
position in a different physical state (¢.g., crystalline and 
glassy )—this is experimentally feasible only in the more easily 
crystallized batches. 

4. Lack of further change when held at the same tempera- 
ture for a longer time—this is the best criterion. 

5. Uniformity of distribution of the phases—this is useful 
to detect inhomogeneity in the original sample. 


None of these criteria is really adequate to detect differences 
between stable and metastable equilibrium when the change is 
sluggish as in the extension of the field of tridymite over that 
of quartz, and the fourth is of no value in the high alkali 
preparations where the composition is changing during the 
duration of the run. The third criterion was applied in many 
cases but only the data for the crystalline material are given 
in the following tables to save space. Initiation of crystalliza- 
tion was so sluggish in most cases that nuclei (traces of 
crystalline material) were added to all quench runs starting 
with glassy material. As with many silicate systems the de- 
termination of the appearance of a phase at equilibrium 
practically requires that that phase be present, at least in 
traces, in the starting material, whether it be mainly crystal- 
line or glassy. For this reason compositions having compli- 
cated crystallization paths involving the formation and sub- 
sequent disappearance of several phases during the cooling 
history were crystallized by various methods on separate por- 
tions, so that starting material having all the necessary 
phases present to establish a given phase change was available 
for each such determination. Obvious cases of minor local 
inequilibrium such as failure of a crystal to react with liquid 
owing to its inclosure by another crystal have been ignored 
in the data tables. Approximately 1300 quenching runs were 
made at temperatures which were not immediately adjacent 
to the temperature of a phase change and hence are not listed 
in the tables of quenching data. These served, however, as con- 
firmatory evidence to substantiate those runs made at the more 
critical temperatures. In those cases in the data tables where 
two runs establishing a phase change have been made for dif- 
ferent lengths of time, confirmatory runs at nearby tempera- 
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tures during the course of the work had established the shorter 
length of time as adequate to obtain equilibrium. 

Temperature Measurement.—Temperatures were determined 
with thermocouples of Pt against Pt 90 Rh 10 measured with 
a shielded, double White potentiometer and galvanometer 
system reading directly to microvolts. The annealed thermo- 
couples (Roeser and Wensel, 1935) were calibrated as fre- 
quently as drifting necessitated.** The melting points of pure 
NaCl (800.4°), gold (1062.6°), diopside (1891.5°), and pseu- 
dowollastonite (1544°) were used to establish the temperature 
scale. The NaCl was prepared for use by two precipitations 
with concentrated HCl from a solution of high purity reagent 
grade NaCl, followed by washing, drying, grinding, and final 
drying at 100°. It had a negligible melting range. The gold 
was “proof gold,” obtained from the U. S. Mint, and also had 
a negligible melting range. This was not true of the diopside and 
pseudowollastonite, which, though synthesized from the purest 
available chemicals, had appreciable melting ranges (about 
30 and 10 microvolts respectively from the first visible liquid 
formation to the disappearance of the last crystal). For this 
reason the temperature of melting of the last trace of crystals 
was considered to be the melting point. These four melting 
points are thus the basis of all of the temperature measure- 
ments, the deviation of the thermocouple from standard ref- 
erence tables being used for convenience in interpolation 
(Sosman, 1910). 

Many of the possible sources of error in the correct measure- 
ment of temperatures are eliminated or minimized by perform- 
ing the calibration runs in exactly the same manner as the 
quenching runs,’® so that thermal gradients within the sample 

18 The drift between calibration runs is variable and only its order of 
magnitude is predictable. High alkali compositions and high temperatures 
considerably increase the rate of drift. A platinum quench envelope sticking 
to a thermocouple necessitates immediate recalibration. The use of iridium- 


free platinum throughout the apparatus aids a great deal in stabilizing 
the emf of the thermocouples. (See Day, Shepherd, and Wright, 1906, p. 
286.) 

19 In the case of gold calibration, a tiny sliver of gold sheet is held in 
a thin porcelain tube 5 mm. long with asbestos plugs. This tube is then 
inserted in a platinum envelope similar to those used in the quenching 
’ runs. After the calibration run the sliver is examined; if the temperature 

was below the melting point of gold, the sliver will be unaffected; if above, 
it will have coalesced to a globule. The coalescence from a white powder 
to a clear mass indicated melting in the case of NaCl. The diopside and 


pseudowollastonite runs were examined microscopically for amounts of 
glass and crystals. 
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and thermocouple are identical in both cases. The effect of in- 
homogeneities in the thermocouple wire away from the junction 
itself and differences in the temperature of the thermo- 
couple junction and the sample itself thus become relatively 
insignificant. All cold junctions were kept in melting ice. 

The hottest portion of the furnace varies with temperature 
and time and hence must be located prior to the run by ex- 
ploring the furnace tube with the thermocouple; the thermal 
gradients and sample size are such that no portion of the 
sample may vary more than a few tenths of a degree from 
that of the thermocouple junction, the method of calibration 
correcting for a portion of this. The details of furnace con- 
struction and regulation have been published previously 
(Adams, 1924; Roberts, 1925, 1941; and Faust, 1936). Fur- 
nace temperature control can be considered to be constant to 
+ 0.8° approximately for runs of a few hours duration and 
+ 0.5° approximately for longer runs. The temperature dur- 
ing runs of over six hours at temperatures of over 1400° was 
only constant to + 1.0°.” 

The frequency of calibration of the thermocouples, and hence 
the accuracy to which any given temperature determination 
on the following pages was made, varied with the sensitivity of 
the phase change being measured,” the temperature, and the 
nature of the batch (change in composition during the run 
makes accurate temperature measurement superfluous in the 
high alkali batches). Relative temperature determinations 
made by running several samples simultaneously in the same 
furnace are, of course, far more accurate than absolute deter- 
minations. By this means small differences between batches 
could be detected although they were less than the absolute 
accuracy of the temperature measurements as determined by 
the calibration technique. The frequency and nature of the 
calibration and the spacing of the quenching runs were such 
that most of the liquidus determinations are accurate to + 3°. 
Whenever possible, more careful work was done in the vicinity 

20 This is due mainly to volatilization of platinum from the winding of 
the furnace and consequent change in its electrical resistance, since the 


furnace temperature control is based on the furnace winding being one 
“arm” of a balanced Wheatstone bridge. 

21 Theoretically phase changes should occur with an extremely small 
AT, but in silicate systems of this type the rate of attainment of equilibrium 
makes such determinations experimentally impractical when AT approaches 


1° or less, except for the more fluid batches and at the melting points 
of compounds. 
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of the invariant points, with the resultant accuracies as stated 
in table 5. 


LIMITING BINARY SYSTEMS 


Three binary systems limit a ternary system. Two of the 
limiting binary systems, MgO-SiO, (Kanolt, 1913; Ruff, 
1913; Bowen and Andersen, 1914; and Greig, 1927), and 
K,0-SiO, (Niggli, 1913; Morey, 1914; Morey and Fenner, 
1917 ; Kracek, Bowen, and Morey, 1929, 1987 ; Kracek, 1930a, 
1932; Goranson and Kracek, 19382; and Kréger and Fingas, 
1933), have been studied and the diagrams of these are re- 
produced in figures 1 and 2. All of the compositions prepared 
in this study were ternary mixtures, hence no new data was ob- 
tained on these limiting binaries. Extrapolation from the 
ternary data to the sidelines, however, is concordant with the 
data on these two binary systems. No data on the system 
K,0-Mg0O have been published, and since all compositions made 
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Fig. 1. Equilibrium diagram of the system MgO-SiO,, after Bowen and 
Andersen (1914) as modified by Greig (1927). 
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in this investigation had over 30 per cent SiO, in them, nothing 
can be said here concerning it. 


EXPERIMENTAL RESULTS 


General Statement.—The compositions by synthesis of the 
133 mixtures studied are listed in table 1 and plotted on figure 
8, along with the boundary curves (heavy, solid or dashed 
lines), fields (primary phase areas), and tie lines (light lines). 
Arrows on boundary curves indicate direction of falling 
temperature. Figure 4 is a plot of the isofracts and curves of 
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Fig. 2. Equilibrium diagram of a portion of the system K,O-SiO,, after 
Kracek, Bowen, and Morey (1937). 
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equal density for the various glasses obtained. The data for 
the indices of refraction are given in table 1. No great ac- 
curacy is claimed for these data; they were obtained by the 
immersion method, with white light, on the unannealed glasses 
quenched from temperatures above the liquidus. They were use- 
ful, however, in determining very roughly the composition of 
the liquid in some quench runs and as convenient checking de- 
tails. Most of these indices are probably accurate to + 0.003, 
although a few, particularly in the higher alkali portion where 
water adsorption and reaction with the index liquid causes 
trouble, may have less accuracy. Reaction with the index oil 
prevented obtaining a figure in several cases. The sharp bend 
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Fig. 8. Equilibrium diagram of the system K,O-MgO-SiO,, showing the 
compositions studied (plotted in weight per cent, by synthesis), compounds, 


boundary lines, and composition triangles (Alkemade lines). The tridymite 
field below 867° is metastable. 
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in the isofracts and in the curves of equal density” at the 
1:1 K,0:Mg0O ratio line of compounds is particularly inter- 


WEIGHT PERCENT 


Fig. 4. Index of refraction and density of the glasses. See text for de- 
tails. Circles indicate the composition of compounds. 


22 These curves of equal density are merely approximate as only 12 
glass densities were determined; these data, although inaccurate for 
several reasons, were adequate to establish the curvature shown. The actual 
densities found were as follows: 
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esting in view of recent work on the effect of the molecular 
structure of glass on its properties (Lukesh, 1948), although 
more accurate refraction determinations as well as density 
measurements are needed. 

In figure 5 the isotherms and temperatures of invariant 
points have been added. The essential data from which this dia- 
gram has been drawn are given in tables 1-5. Compositions 
higher in MgO or K,0 than those in the tables were not studied 
because of the lack of direct geologic interest and because of 
experimental difficulties. All data below the join 2MgO-Si0,- 
K,0°2Si0, (LE, fig. 5) are considered to be of a reconnais- 


76 
WEIGHT PER CENT 

Fig. 5. Equilibrium diagram of the system K,O-MgO-SiO,, showing 
isotherms and invariant points. The tridymite field below 867° is meta- 
stable, hence binary eutectic point B and ternary eutectic point R are 
both metastable. For an enlarged plot of the sub-system Forsterite- 
K,O.Mg0O.5SiO,-SiO, (triangle LWS’) see figure 8. Data below a line 


between 2MgO-SiO, and K,0.2SiO, (LE) are of a reconnaissance nature; 
see footnote 49. 
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sance nature, although they are fairly concordant among them- 
selves. In addition to the compounds already known to occur 
in the limiting binary systems MgO-SiO, and K,O-SiO,, four 
new ternary compounds were found during the course of the 
work. They have the compositions K,0-Mg0-SiO,(?), K,0- 
Mg0-3Si0., and K,0°5Mg0'12Si0,. A por- 
tion of figure 5 (the triangle LWS’, 2MgO-Si0.-K,0:MgO- 
5Si0,-SiO,) has been recalculated in terms of these three com- 
ponents and replotted on a larger scale in figure 8 to show 
courses of crystallization, etc. 


2000 


|2MgO-Si0, 


WEIGHT PER CENT 


Fig. 6. Sections through the ternary system along the various binary 
sub-systems. The two cross-hatched lines indicate non-binary (ternary) 
equilibrium for those portions. Several of the binary systems that are almost 
wholly non-binary at liquidus temperatures have been omitted. Data taken 
from tables 1 to 3. 
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COMPOSITION OF THE PREPARATIONS STUDIED 


The percentages are weight per cent by synthesis (see text and 

table 7 for analysis of high alkali batches). The molecular weights 

used in the calculation were K,O-94.20, MgO-40.32, and SiO,-60.06. 
Mol ratio is the molecular ratio K,0:MgO:SiO,,. 


Mol Per Cent nof Mol Per Cent 
Ratio Prep. K,O0 MgO SiO, Glass Ratio Prep. 


1:1:6 
1:2:6 
13:6 
1:4:6 
1:5:6 
1:6:6 
1:8:6 


19.03 8.15 7282 1.499 
17.60 15.07 67.38 
1636 21.01 62.63 
15.29 26.19 58.52 
1435 30.72 54.93 
13.53 34.73 51.74 
12.12 41.51 46.37 
1:12 36.99 15.83 47.18 
1:22 $1.94 27.34 40.72 
123 26.58 22.71 50.76 
1:1:8(2) 29.938 1281 57.26 
18.24 12.00 69.76 
$7.36 8.60 54.04 
26.57 5.68 67.75 
25.14 10.76 64.10 
18.65 10.00 7135 
8.00 23.00 69,00 
$3.09 14.16 52.75 
7.00 18.00 75.00 
13.00 10.00 17.00 
40.17 860 51.23 
26.00 14.00 
2787 641 

$5.00 5.00 
21.67 9.27 

22.69 19.43 

10.00 10.00 

7.00 15.00 

11.00 13.00 

19.00 

41.76 

34.00 

29.49 


2 


7.00 
14.38 
16.00 
11.00 

7.00 

4.00 
15.00 
22.00 
21.00 
24.00 
10.00 
18.00 


98 
nof 
1 SiO, Glass 
: | 12.00 9.00 79.00 1.492 
16.00 12.00 72.00 1.499 
18.00 16.00 71.00 1.507 
2.00 26.00 72.60 1.534 
26.00 2.00 72.00 1.4938 
27.18 2.00 170.82 1.496 
31.00 2.00 67.00 1.497 
28.46 15.00 56.54 1.515 
2:2:9 23.27 9.96 66.77 1.500 
4.00 31.00 65.00 1.551 
$1.68 16.00 52.32 1.517 
40.12 5.00 54.88 1.508 
$2.00 6.00 62.00 1.508 
$1.50 4,00 64.50 1.499 
27.50 8.50 64.00 1.500 
25.00 9.50 65.50 1.501 
24.48 7.50 68.02 1.497 
22.50 9.50 68.00 1.498 
28.99 7.50 68.51 1.497 
20.02 10.50 69.48 1.498 
19.50 6.00 74.50 1.493 
i 17.00 6.50 76.50 1.490 
; 22.91 17.50 69.59 1.497 
* 42:17 170 2549 545 69.06 1.498 
4 1:2:5 71 19.83 1697 63.20 1.515 
1:29 72 18.17 11.27 75.56 1.498 
73 11.00 8.50 80.50 1.490 
: 14 15.00 7.00 78.00 1.490 
i 15 24.00 2.50 73.50 1.493 
t 76 28.61 3.20 68.19 1.497 
ry 17 22.53 1146 66.00 1.503 
78 21.00 12,00 67.00 1.503 
79 19.72 12.05 6823 1.502 
80 5.50 12.50 82.00 1.496 
11.00 82.00 1.494 81 9.00 8.50 82.50 1.488 
1:32:38 1230 73.32 1.504 82 16.00 6.50 77.50 1.490 
8.00 76.00 1.493 83 18.50 5.50 176.00 1.490 
5.00 84.00 1.4838 84 24.50 3.50 172.00 1.496 
9.00 84.00 1.487 85-85a 25.71 3.50 70.79 1.496 
17.00 79.00 1.505 86 27.00 3.50 69.50 1.497 
5.00 80.00 1.485 87 17.50 9.00 738.50 1.497 
400 174.00 1.495 88 16.50 11.00 72.50 1.498 
8.00 71.00 1.497 89 17,00 11.50 71.50 1.498 
5.00 71.00 1.497 13:10 90 11.55 14.88 73.62 1.608 
7.00 83.00 1.485 91 837.78 5.64 56.78 1.508 
7.00 80.00 1.487 92 84.14 8.97 5689 1.508 
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Mol Per Cent 
Ratio Prep. K.0 MgO 


3:2:6 


2:13 


1:1:1 121-12la 

2:12 122 
128 
124 


110 00 51.19 1.515 
111 


‘ .00 40.00 1,531 
112 J 38.98 1.528 
113 72.50 1.509 


(1) These batches would not melt completely within the limit of the 
furnaces (1600°). 


(2) This batch has the composition of the compound K,O.Mg0O-3Si0,. 

(3) This batch has the composition of the compound K,0.M .5SIO,. 

(4) This batch has the composition of the compound K,0.5Mg0O. 12Si0,. 

(5) This batch has the composition of the compound K,O. MgO. Si0,. 

As may be seen in table 1, a number of the compositions 
were made to have integral molecular ratios of the component 
oxides. This procedure helped considerably to locate the com- 
position of the compounds during the early part of the work 
in that it eliminated many possible combinations from further 
consideration. Other compositions not having simple unit 
weight per cent ratios were made to fall on the various binary 
joins within the ternary system. Of these joins, those which 
illustrate partial or complete binary equilibrium at liquidus 
temperatures have been plotted in figure 6. 

Whenever possible boundary-curve determinations were made 
(table 3), although the optical properties and crystallization 


habit of some of the yhoo precluded experimental determin- 
ations in certain cases.” 


23 For example, the boundary curve 4’B’ (fig. 5) could be determined 
only on batches lying in the K,O.MgO-5Si0, field as the tiny isotropic 
crystals of this phase in a glass of almost the same index of refraction 
were impossible to see if much of the finely divided, higher index, bi- 
refringent K,O.MgO-3SiO, was. present. 


Mol Per Cent nof nof 

Ratio Prep. K.0 MgO SiO, Glass SiO, Glass 
938 26.00 7.50 66.50 1.500 114 $0.00 2.00 68.00 1.497 
94 24.00 13.50 62.50 1.508 115 30.00 4.00 66.00 1.499 
95 23.00 13.00 64.00 1506 [i= 116 $9.06 11.14 49.80 1.516 
96 19.84 13.31 6685 1.505 117 41.54 5.50 52.96 1.510 
97 20.45 11.00 68.55 1.500 118 48.03 18.00 48.97 1.522 

1:5:12 98-98a(4) 9.27 19.83 70.90 1517 Him 119 46.08 986 44.06 1.522 
99 22.08 945 6847 1.496 120 45.80 25.00 29.20 (1) ts 
100 22.90 9.81 67.29 1.499 48.41 20.72 80.87 1.547(5) 
101 23.73 10.16 66.11 1.500 54.01 1156 3448 1.527 nt 
102 24.56 10.51 6493 1.501 59.23 8.00 37.77 reacts 
108 26.21 11.22 62.57 1.504 87.80 14.00 48.20 1.519 
104 $0.00 5.81 64.19 1.501 125 41.10 650 5240 1.512 
105 28.50 12.00 64.50 1.504 126 $6.39 11.63 51.98 1.514 
106 14.50 8.00 17.50 1.490 127 35.00 14.00 51.00 1.517 
107 21.50 56.00 78.50 1.492 128 46.78 138.12 40.10 1.525 
108 13.00 8.00 79.00 1.490 129 60.46 1.00 38.54 reacts 
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DATA ON LIQUIDUS DETERMINATIONS 


Abbreviations Used: L-large; M-moderate; S-small; VS- 
very small; R-rare; VR-very rare; VVR-exceedingly 
rare—right at the temperature of the phase change. These 
are rough visual estimates under the microscope, and are 
useful for comparison purposes only, as crystal habit 
and index of refraction may affect this estimate seriously. 


Above Liquidus Below Liquidus 
Quenched sample Quenched sample contained 
all glass crystals in glass 
Prepara- Amount of 
tion Temp. °C Time Temp. °C Time crystals 


Primary crystal phase—Cristobalite, SiO, 


1535 5 hrs. 1530 4 hrs. — R 
1590 3” 1580 


Primary crystal phase—Tridymite, SiO, 


930 3 days 
1325 20 hrs. 
1485 
1115 days 
1320 
1135 2 
1210 . hrs, 
985 days 
1340 hrs. 
785 days 
760 (2) 


Primary crystal phase—“Clinoenstatite,” MgO-SiO, (3) 


1480 1475 5 
1310 1305 4 
1465 1460 6 
1360 1355 17 
1420 1415 17 
1265 1260 4 
1245 1240 22 
1380 1375 19 


Primary crystal phase—Forsterite, 2MgO-SiO, 


1805 5 hrs. 1300 5 hrs. 
1465 (4) ” 1460 

1545 1540 
1590 
1210 
1475 
1125 
1490 
1180 
1445 
1880 
1165 
1530 


< 


<< 


” 


< 


= 
eanBancen 


100 
40 
30 R 
38 s 
38 VR 
41 VR 
45 Ss 
46 Ss 
73 R 
4 
81 VR 
i 110 VR(1) 
ill VR(1) 
5 19 hrs. 
27 days 
28 hrs. 
| 29 
35 
36 days 
72 hrs. 
90 
2 
3 
5 
8 
10 
12 8 
17 
22 VV 
26 
49 
54 
56 
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Taste 2 (Continued) 


Above Liquidus Below Liquidus 
Quenched sample Quenched sample contained 
all glass crystals in glass 
Prepara- 


tion Temp. °C Time Temp. °C 


57 1200 5% hrs. 1195 

71 1355 days 

78 1065 

79 1095 

94 1125 

95 1110 
1190 
1470 
1465 
1540 
1410 
1365 
1395 


Primary crystal phase—K,O0-5M 
1175 1 day 1170 
1020 44,” 1015 
1165 1160 
1135 49 1130 
1140 1185 
1055 = 1050 


1100 hrs. 1095 12 hrs.(7) 
crystal 
4 days 


” 


days 


brs. 
” 
»” 
» 


bo 


Sm mm 


wer 


1060 


Dm 


Ly 
7 
5 
5 
1 
6 
2 
5 
2 
8 
8 
3 
2 
3 
4 
4 
5 
2 
5 
7 
6 
6 
5 
5 
4 
5 
3 
5 


ae nore 


101 
P 
of 
crystals 
hrs. R 
R 
days VR 
R 
days VR 
VVR 
hrs. R 
(5 R 
(5 
(6 R 
R 1 
» q 
R 
-12Si0, 
20 days 
37 Ss 
47 VR 
88 
89 VR 
106 
108 
1 
14 980 975 
23 980 975 
25 1090 1088 ”(8) Fi 
42 885 880 ” i 
43 1080 1075 
44 970 965 4 
55 1080 1075 “ng 
60 790 785 
62 1055 1050 i 
63 1055 1050 
64 1085 1084 7 
65 1060 1055 x 
66 1070 1065 bs 
67 1010 1005 53 
68 1005 1000 be 
69 1065 1060 iy 
70 985 980 
75 765 760 
76 835 830 
77 1050 1045 
82 995 990 ? 
83 975 970 af 
84 860 855 * 
85 865 860 by 
85a 865 860 . 
86 860 855 na 
87 1045 1040 - 
93 1030 1025 ge 
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Taste 2 (Continued) 


Above Liquidus 
Quenched sample 
all glass 
Prepara- 
tion Temp. °C 


1065 
1084 
1080 
1070 
1050 
1035 

955 
1005 
1075 
1045 


38 


eww @ 

ee 
3 < 


crystal 
day 1138.5 2 days(10) 
hrs. 


days 


hrs. 1050 
days 1030 
1130 
880 

1090 

930 

995 

960 

1078 

1075 

880 


1000 
965 
1080 
1080 
885 


Primary crystal phase—K,Q-MgO.-SiO,(?) 


1090 days 1088 5 days 
1270 hr. 1265 1 hr. 
1385 min. 1380 15 
1385 1380 15 
1470 15 
1220 hr. 1 
1310 min. 15 

985 hr. 10 
1140 day 1 
1210 hr. 1200 1 
1450 min. 1440 15 

905 hr. 900 1 


< 


Primary crystal. phase—Periclase, MgO 
hr. 1595 hr. 


Primary crystal phase—K,O.4Si0, 
12 days 735 7 days 
9 


2 ” 
5 ” 


102 
Below Li 
crystals in glass 
Amount of 
crystals 
97 
99 
100 
101 
102 
105 
107 
109 
112 
118 
115 845 
Primary 
1135 
13 1055 
15 1035 5 
18 1133.5 2 
24 885 be 
83 1095 20° hrs. 
59 935 5 days 
61 
91 20 hrs. 
92 2 days 
108 
21 
116 
118 
119 
122 
123 
124 
; 125 
126 
127 
128 
130 
52 739 
75 725 720 §(11) 
mS 
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Taste 2 (Continued) 


Above Liquidus Below 
Quenched sample Quenched sample contained 
all glass crystals in glass 
Prepara- Amount of 
tion Temp. °C Time Temp. °C Time crystals 


Primary crystal phase—K,O.2Si0, 


31 995 
860 
739 


Primary crystal phase—K,O-SiO, 
980 10 hrs. . 970 1 hr. 


(1) Metastable, only tridymite in the starting material. 

(2) Held 4 days at 740° first. 

(3) May include other modifications—see discussion in text. 

(4) Quenched glass contained exceedingly rare crystals; right at the 
liquidus temperature. 

(5) These preparations have the composition K,O.5MgO.- 12SiO,. 


(6) This preparation was originally of the composition K,O.MgO-.Si0, 
—see data on alkali losses. 

(7) This result was duplicated in a run at 1090° for 4 days. 

(8) This preparation has the composition K,0-Mg0.5Si0,. 

(9) “Hydrothermally” crystallized material used. “Dry” crystallization 
yielded only K,O0-4SiO, crystals and hence a metastable liquidus 
determination. 


(10) This preparation has the composition K,O.MgO-SSiO,. 
(11) Metastable; see footnote (9) above. 


Stability Fields for the Various Crystalline Phases of the 
Limiting Binary Systems.—The fields of stability of the prim- 
ary crystalline phases within the ternary system are shown in 
figures 3 and 5 and are labeled in figure 3. 


Cristobalite——The field of cristobalite (IS’ANI, fig. 5) ex- 
tends from 1470° up to the melting point of cristobalite at 
1713 + 5° (Greig, 1927, p. 12). The crystals obtained at 
higher temperatures were usually rounded, though occasionally 
facetted masses showing very low birefringence and the char- 
acteristic complex twinning caused by the 8 > < inversion 
were formed. The smaller crysals formed at lower temperatures 
appeared sensibly isotropic because of their size. 

In that portion of the cristobalite field lying between H’H 
and the dashed curve in figure 5 two conjugate immiscible 
liquids are in equilibrium with cristobalite at appropriate 
temperatures between the temperature of H’H and that of the 


108 
hrs. 990 hrs. 
5 days 855 4 days R 
58 945 940 17 hrs. 
117 930 l4y,” 925 8 
131 920(4) hr. 910 hr, 
| 
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DATA ON BOUNDARY CURVES 


Abbreviations used: 
Phases—Cris Cristobalite, SiO, Amounts— 


Trid ‘Tridymite, SiO, L 

Clino Clinoenstatite, MgO-SiO,(1) M 

Forst Forsterite, 2MgO-.SiO, Ss 

Peric Periclase, MgO very small 

K,O-MgO-SiO, rare 

1:13 K,0-Mg0.3Si0, very rare 

K,O-Mg0-5SiO, exceedingly rare, 
1:5:12 K,O-5Mg0.12Si0, right at 

KS, K,0.-4Si0, temperature of 


KS, K,O.2Si0, the phase change 
KS K,0.-SiO, 


gl glass 


Preparation Temperature Time Phases present in quenched 
sample (2) 


Between Cristobalite and Clinoenstatite 


1510 6 hrs. M cris 

1505 - M cris and VR clino 
1535 = M cris 

1530 oy M cris and § clino 
1475 All gl 

1470 43 S cris and R clino 


Between Tridymite and Clinoenstatite 
1240 4 days L clino 
1235 axa L clino and R trid 
1400 5 hrs. M clino 
1395 i L clino and R cris (metastable) 
1380 L trid 
1375 2" L trid and R clino 
1175 7 days M trid 
1170 , M trid and VVR clino 
1275 Blin M trid 
1270 G:* M trid and S clino 


Between Tridymite and K,0-5Mg0O- 12SiO, 
1140 2 days L trid 
1135 L trid and § 1:5:12 
1050 L trid 
1045 a L trid and VR 1:5:12 
1135 L 1:5:12 
1130 L 1:5:12 and R trid 
975 L 1:5:12 
970 L 1:5:12 and R trid 
1075 L 1:5:12 
1070 L 1:5:12 and VR trid 


Between Tridymite and K,O-MgO-5SiO, 
965 L trid 
960 L trid and R 1:1:5 
93a L 1:13 
925 L 1:1:5 and R trid 


Norenwe w 


104 
40 
50 
80 
27 
35 
; 39 
4 73 
81 

i 45 
46 

47 

106 

108 

41 

67 


The System K,0-Mg0-SiO, 


Taste 8 (Continued) 


Preparation Temperature Time Phases present in quenched 
sample (2) 


21:5 
:1:5 and VR trid 
id 


trid and M 1:1:5 
1:1:5 

1:1:5 and M trid 
1:1:5 

1:1:5 and § trid 


Between Tridymite (metastable) and K,0.4Si0, 


111 750 4days R trid 
745 4” R trid and R KS, 


Between Clinoenstatite and K,0-5MgO-12SiO, (Reaction curve) 


17 1175(4) 7 days L clino and VR 1:5:12 
1172(4) . M clino and M 1:5:12 
29 1170 L clino 
1165 L clino and VR 1:5:12 
1160 L 1:5:12. 
36 1175 L clino 
1170 R clino and L 1:5:12 
1165 L 1:65:12 
1165 L clino 
1160 L clino and M 1:5:12 
1175(4) L clino 
1172(4) R clino and L 1:5:12 
1170(4) L 1:5:12 
1175 (4) L clino 


Starting 
material 
1170(4) Lclino and VR 1:5:12}“bombed” 
1:5:12 
crystals 


Between Forsterite and Clinoenstatite (Reaction curve) 


L forst and VVR clino 
L forst and S clino 

L clino and R forst 

L clino 

All glass 

S clino and R forst 
M clino and R forst 
forst 

L forst and L clino 

L clino and VR forst 
L clino 

S forst 

M forst and § clino 
VR forst and L clino 
L clino 

S forst 

M forst and R clino 
L clino and VR forst 
L clino 


Dp 


ss3ss3 3 3 


ee 


105 
68 980 5 days L 1 
74 980 Rt 
975 5” M 
82 990(3) 6 ” R 
985(3) 5 ” M 
83 925 16 hrs. L Hig 
920 4days L 
| 
17 1460 6 hrs. 
1455 6” 
1395 26” 
1390 4 days 
48 1180 
1175 
1160 
49 1285 
1280 
1200 
1195 
56 1530 : 
1525 
1485 
1480 
98 1455 
1450 
1400 
1395 


106 Edwin Woods’ Roedder 


Tastz 8 (Continued) 


Preparation Temperature Time Phases present in quenched 
sample (2) 


98a 1445 M forst 
1440 M forst and M clino 
1400 L clino and VVR forst 
1890 L clino 
134 1890 R forst 
1385 S clino 
Between Forsterite and K,O-5MgO-12SiO, (Reaction curve) 
2 1100 L forst 
1095 L forst and R 1:5:12 
12 1120 R forst 
1115 VR forst and S 1:5:12 
1110 VR forst and M 1:5:12 
1105 M 1:35:12 
1120 L 1:5:12 and § forst 
1115 L 1:5:12 
1070 M forst 
1065 S forst and M 1:5:12 
1057 L 1:5:12 


3 
& 


3338 


aaa 


Between Forsterite and K,O-Mg0O-5Si0O, 
Sdays L 1:18 

L 1:1:5 and M forst 

R forst 

M forst and M 1:1:5 

M forst and L glass 


L forst, L 1:1:5, and no gl 
1:13 
L 1:1:5 and S forst 
L 1:1:5 and L glass 
1055(5) L 1:1:5, R forst, and L gl 

1050(5) L 1:1:5, L forst, and VVR gl 
Between Forsterite and K,O-Mg0O-3Si0O, 

1125 19 hrs. L forst 

1120 2 days L forst and M 1:1:3 

1115 L 1:13 

1110 L 1:1:3 and L forst 

1100 M forst 

1095 M forst and R 1:1:3 

1182(6) 2 days L forst 

1130(6) fb L forst and L 1:13 

1125 19 hrs. L forst 

1120 2 days L forst and M 1:1:3 

1050 L forst 

1045 ies L forst and § 1:1:3 


Between Forsterite and K,O.MgO-SiO, (6a) 
1110 4 days L forst 
1100 L forst and L 1:1:1 
1260 L forst 
1255 > L forst and L 1:1:1 


Between Periclase and Forsterite (Reaction curve) 
1525 2 hrs. M peric and VR forst 
1490 a -R peric and R forst 
1480 M forst 


1057 (5) 


Ssss3333 8 


| 
| 
78 
96 
105 
112 | 
i 
10 
} 18 
2 
f 
57 
94 
8 
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The System K,0-Mg0-Si0, 
Taste 8 (Continued) 


Preparation Temperature Time Phases present in quenched 
sample (2) 


Between K,0-.MgO.-SiO, and K,0- SiO, (6a) 
920 1 hr. L 1:1:1 


122 
910 2 days L 1:1:1 and R KS 
128 935 1 hr. L 1:1:1 


925 4 hrs. L 1:1:1 and § KS 


Between K,O.MgO-SiO, and K,0-2Si0, (6a) 
910(7) L and L gl 


2 
905 (7) L 1:1:1, S KS,, and gi 
910 2 L 
1 
1 
1 
1 


905 L 1:1:1 and L KS, 
730 L 


720 L 1:1:1 and S KS, 
8 K 
905(7) 1 L KS, and R 1:1:1 


Between K,O-MgO-SiO, and (6a) 

1:1: and L 1:13 
: 


. 


K,0.2Si0O, 
1:33 

:1:3 and S KS, 
1:3 

1:3 and R KS, 


915 
910 
875 M 
870 L 
920 L 
915 L 
915 L 
910 L 
L 
L 
L 
L 
L 


L 
L 
L 
L 
L 
L 
L 
L 
an 
L 
L 


1 
1 
1 
1 
1 
1 
1 
d 
1 
1: 


1 
1: 
KS, 
KS, and R 1:1:3 
KS, 
KS, and R 1:13 
820 
815 1: 
985 (8) 1 
931(8) 
935 (8) 
931(8) 


Between K,O MgO ssiO, and 
days 


1:3 and R KS, 

1:33 

1:1:3, L KS,, and R gl 
1:1:3 and M gl 
1:1:8, L KS,, and R gl 
K, 


O.MgO.5Si0, 


4 
5 
3 
3 
3 
0 


» 
” 
» 
” 
» 
” 
» 
» 
» 
” 


5 
4 
6 
0 
5 
6 


119 
125 
130 
131 
116 
118 
124 
1120 and R 1:1:3 
127 1140 
1185 1:1 and L 1:13 
Between K,O.MgO-3Si0, 
18 
24 j 
$1 
58 
59 
91 
92 
55 1:5 and L gl 
1030(9) 1:5, M 1:1:3, and R gl 
62 135 
1025(9) 1:5 and § 1:13 
64 1025 1:5 and R gl 
1020 1 1:5, S 1:1:3, and R gl 
98 860 1:5 
855 1:5 and R 1:13 
100 1025 1:5 and S gl 
1020 1 1:5, L 1:1:3, and R gi 
101 1080(9) and L gl 
1025 (9) 1:5, S 1:1:3, and L gl 
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Taste 8 (Continued) 


Preparation Temperature Time Phases present in quenched 
sample (2) 


102 1030(9) 5 days 
1025(9) 6 ” ‘5, R 1:1:3, and L gl 
118 1025 21:5 
1020 21:5 and L 1:1:3 


Between K,O.MgO.5Si0, 


810 
805 
785 
780 
745 
740 1:1:5 and R KS, 


Between K,O-.2SiO, and K,O.4SiO, 


710 5 days S KS, 
705 S KS, and R KS, 
695 = L KS, 
690 S..* L KS, and § KS, 


Between K,O-MgO-5SiO, and K,O0-5MgO. 12SiO, 


L 1:5:12 
L 1:5:12 and S 1:1:5 
All glass 
S 1:1:5 and § 1:5:12 


1065 (10) 
1060(10) 

1040 

1085 

1040 

1035 

1050 

1045 

97 1060 

1058 


:1:5 and § 1:5:12 
21:5 

:1:5 and VR 1:5:12 
25:12 

25:12 and VR 1:1:5 


5:12 and M 1:1:5 
21:5 

and § 1:5:12 
(1) May include other modifications—see discussion in text. 

(2) Glass is present in all cases unless otherwise stated. 

(3) These runs determine binary eutectic point Q. 

(4) These runs determine binary reaction point T. 

(5) These runs determine binary eutectic point Y. 

(6) These runs determine binary eutectic point I’. 


(6a) The isotropic crystals obtained in some of these runs may not be the 
compound K,0.Mg0O-SiO,, see footnote 49. 


(7) These runs determine the minimum temperature of the binary eutectic 
point L’. 


(8) These runs determine binary eutectic point E’. 
(9) These runs determine binary eutectic point A’. 
(10) These runs determine binary eutectic point V. 


WEL 


and K,0-2Si0, 
23 
60 
Ss 
115 
53 
114 
12 1060 5 days 
1055 
16 1055 Behe 
1050 
$7 1000 S 1:5:12 
995 S 1:5:12 and 1:1:5 
L 
L 1 
” Ss 1 
” M l 
Ll 
L1 
L 1:5:12 
| 


The System K,0-Mg0O-SiO, 


TaBLe 4 


DATA ON TERNARY INVARIANT POINTS 


(Abbreviations given at heads of tables 2 and 3. The com- 
positions of the ternary points are given in table 5.) 


Prepara- Tempera- 
tion ture °C 


Time 


Phases present in quenched sample (1) 


Reaction Point N—Cristobalite, Tridymite, Clinoenstatite, and liquid 


80 


1475 
1470 


6 hrs. 
6 » 


All glass 
S cris, R clino, and gl 


Reaction Point O—Tridymite, Clinoenstatite, K,0.5MgO-12Si0,, and liquid 


27 


45 


47 


111 


1170 
1165 
1160 
1140 
1135 
1165 
1160 
1135 
1130 
1170 
1165 
1160 


7 days 


L cline, S trid, and gl 

S clino, M trid, L 1:5:12, and gl 
VVR clino, M trid, L 1:5:12, and gl 
L trid 

L trid and § 1:5:12 

All glass 

VR 1:5:12 

L 1:5:12 

M 1:5:12 and R trid 

M trid, VVR clino, and gl 

M trid, M clino, R 1:5:12, and gl 
M trid, S.1:5:12, and gl 


Eutectic Point P—Tridymite, K,0.5MgO-12Si0,, 
K,O.MgO-5SiO,, and liquid 


965 
960 
975 
960 
965 
960 


4 days 
6 ” 


L 1:1:5, R 1:5:12, and M gl 

1:1:5, R 1:5:12, (M trid?), and no gl 
1:1 5, M trid, and L, gl 
1:1:5, L trid, R 1:5:12, and S gl 
1:5:12, R trid, and L gl 
1:5:12, M(?) trid, (Lid :5?), and S-R gl 


L 
M 
L 
M 
M 


Eutectic Point R—Tridymite, K,O.MgO-5SiO,, K,0-.4Si0O,, 
and liquid (metastable) 


165 
760 
725(2) 
720(2) 
720 
710 


7 days 


All glass 

VR 1:1:35 

All glass 

S K,0.4SiO, (metastable) 

L KS,, L trid, and L glass 

L KS, L trid, S(?) 1:1:5, and R gl 


Reaction Point U—Forsterite, Clinoenstatite, K,0.5MgO-12Si0,, and liquid 


48 


49 


88 


1160 
1155 
1150 
1160 
1155 
1150 
1135 
1180 
1140 
1135 


7 days 


» 


R forst, M clino, and glass 

S forst, S clino, S 1:5:12, and gl 

M 1:5:12 and glass 

L clino, M 1:5:12, and gl 

L clino, VR forst, L. 1:5:12, and gl 
R forst, L 1:5:12, and gl 

All glass 

S 1:65:12 

All glass 

VR 1:5:12 


109 
= | 
” 
» 
| 
” 
” 
” 
” 
73 
= 4 
14 5 ” 
4” é 
106 4” id 
4 » 
75 | 
2 
2 ” 
| 
5 ” 
= | 
” 
= 
” 
” 
» 
” 
” 
89 
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Taste 4 (Continued) 


Tempera- 
ture °C Time Phases present in quenched sample (1) 


Eutectic Point X—Forsterite, K,O-5MgO-12Si0,, 
K,O.MgO.5Si0,, and 
1040 4 days L 1:5:12, S forst, L(?) 1:1:5, and L 
1085 4” L 1:5:12, S forst, L(?) 1:1:5, and R gi 
1050 oc S 1:5:12, L 1:1:5, and § 
1045 Pies S 1:35:12, L 1:1:5, R(?) forst, and no gl 
Eutectic Point Z—Forsterite, K,O-MgO-5SiO,, 
K,O-MgO.3SiO,, and liquid 

1015 3 days L forst, L 1:1:5, and L 
1010 L forst, L 1:1:5, L 1:1:8, and VR gl 
1015 L forst, and L gl 
1010 L forst, L 1:1:3, L 1:1:5, and no gl 
1015 M forst, L 1:1:5, and L gl 
1010 M forst, L 1:1:5, L 1:1:3, and R gl 
1015 S forst, L 1:1:5, L 1:1:3, and L 
1010 S forst, L 1:1:5, L 1:1:3, and R gl 

Reaction Point B’—K,O-MgO-5Si0,, 

K,O.-2Si0,, and liquid 

days L 1:1:3, 8 KS,, and M gl 

ad 1:3, S KS,, (M 1:1:5?), and S gl 

8, S KS,, (L 1:1:5?), and R gl 

8, S KS,, (L 1:1:5?), and no gl 


and VR KS, 
5, R 1:1:8, and L gi 
71:5, R 1:1:3, and R KS, 
gl 


810 
805 
800 
790 
785 


ig 
I: 
1: 
780 1: 
I: 
1: 


= 


7 
9 
8 
8 
6 
3 
4 
9 
5 


Eutectic Point K,O.4Si0,, 
K,0-2SiO,, and liquid 
710 5 days 
705 5 
835 6 
830 7 
695 
5 L KS, and S KS, 
4 L KS, and R 1:1:5 
2 L KS,, R KS,, and R 1:1:5 
K,O 


nt K,0.MgO.-Si0,, 
liquid (3) 

, S KS,(?),1:1:1(?), and S gl 
KS, (?), 1:1:1(?), and R gl 
KS,, and M 

S., ($1 and no gl 
Ss 


2h 


48 
48 
48 
42 
48 
2 days 
12 hrs. 
$3 days 
16 hrs. 


$3 days 
12 hrs. 


110 
Prepara- 
79 
97 
17 
95 
105 
118 
59 
60 1:5 
93 
104 885 a 
co 
76 
114 
115 
E 
18 910 
. 905 
$1 910 
905 
58 910 
905 S,, M 1:1:3, (S 1:1:1?), and no gi 
116 910 :1:1, L 1:1:3, and L gl 
900 :1:1, L 1:13, L KS,, and R g 
117 905 KS,, (1:l:1?), and M gl 
900 glass 
124 905 :1:1, L 1:12:38, and R gl 
900 1:1, L 1:1:3, R KS,, and VR gl 
| 


The System K,0-Mg0-Si0, 


Taste 4 (Continued) 


Prepara~ Tempera- 
tion ture °C Time Phases present in quenched sample (1) 


125 905 lhr. L 1I:1:1, L KS,, and L gi 
900 PR L 1:1:1, L KS,, (L 1:1:3?), and R gi 


Eutectic Point G’-—K,0.Mg0.3Si0,, K.0.Mg0-SiO,, 
Forsterite, and liquid (3) 

days L forst and L gl 

L forst, L 1:1:1, and VR gi 

L forst, L 1:1:3, and L g 

L forst, L 1:1:3, (S 1:1:1?), and no gl 

L forst, L 1:1:3, and S gl 

L forst, L 1:1:3, (S 1:1:1?), and R gl 

L forst, L 1:1:3, (S 1:1:1?), and no gl 

L forst, L 1:1:3, and S-R gl 

L forst, L 1:1:3, (S 1:1:1?), and no gl 


(1) Most of these results were obtained by starting with the all. 
material. The intergrowths of very minute crystals with insufficiently 
distinctive optical properties prevented accurate determination of the 
presence of one of the three crystal phases at the ternary point in 
some cases. Reliance then had to be placed on sharp differences in the 
amounts of liquid (glass) formed in samples of the all-crystalline 
material held at various temperatures. 

(2) Only K,O.4SiO, crystals present in the glass at the start. 


(3) The isotropic crystals obtained may not be the compound K,O.MgO. 
SiO,, see footnote 49. 


1110 4 
1100 5 
1110 5 
1105 2 
1110 5 
1105 1 
1095 4 
1110 2 
1105 2 


. 
8 
10 
is 
57 
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SUMMARY OF INVARIANT POINTS 


Point (1) Equilibrium type Crystal phases in equilibrium with liquid at point 


Mel 
elting 


Eutectic Tridymite + K,0-Mg0O.5Si0O, 
Reaction Clinoenstatite + K,O.5MgO- 12SiO, 
Eutectic 


° 
ES 


+ K,O.5MgO- 12Si0, 


» 


+ K. O. 2SiO, 

-3Si0, + Forsterite 
-SiO,(?) + K,O.2Si0, 
-SiO,(?) + K,O.SiO, 
-SiO,(?) + Periclase 
-Si0,(?) + K,O-MgO-3Si0, 


BINARY 


” 


? 
Eutectic? 


Inversion Clinoenstatite +. Cristobalite 4+. Tridymite 
Reaction j K,0O.5MgO. 12SiO, 
Eutectic 

” (3) 
Reaction 
Eutectic 


Reaction K,O. MgO.3Si0, + 
Eutectic K,O-Mg0O.5Si0, + K. 
+ K.O. MgO -Si0, (?) + K,0-2Si0, 
MgO.3Si0, + K,O.MgO.SiO,(?) + Forsterite 
Reaction MgO.SiO, (?) + Forsterite + Periclase 
Eutectic MgO .-SiO,(?) + K,0-2Si0, + K,0.SiO, 


TERNARY 


(1) See figure 5. 


(2) The compositions of the invariant points were obtained mainly by ex- 
tension of the established boundary curves to intersection. The composi- 
tions are in weight per cent, as are the approximate limits. In those 
cases where one component is present in considerably smaller quantities 
than the others, the limit of error is probably less than the figure 
given for that component. 


(3) Metastable. 


112 | 
Ww K,O.MgO-5SiO 
K,0.Mg0.3Si0, 
oO K,O- MgO. SiO, (?) 
Q 
T 
Vv 
Y + Forsterite 
A’ + K,0-MgO.3SiO, 
Db = gO 
E’ gO 
V gO 
L’ ” gO 
N’ = gO 
Pp’ gO 
£0 | 
N 
P id, 
R 
U 
x iO, 
j Z ” 
& B’ 
F’ 
K’ 
M 
| 
| 


The System K,0-Mg0O-Si0, 


Taste 5 (Continued) 
SUMMARY OF INVARIANT POINTS 


Composition (2) 


Approx. 
Point (1) K.O MgO SiO, limits ++ Temperature 
per cent per cent per cent per cent °c 


21.66 9.27 69.06 1089 + 1 
29.93 12.81 57.25 1134 +1 
48,41 20.72 (1650 +) 


15.5 987 + 5 
13.5 1174 +2 
20.0 . 1063 + 3 
20.5 1053 + 2 
24.8 \ . 1030 + 5 
26.5 730 + 5 
38.3 933 + 5 
29.5 1131 +1 
44.8 910 + ? 
60.1 ‘ 950 + ? 

>1650 + ? 

33.0 1110 + 20? 


5.7 1470 + ? 
11.6 J 1165 + 2 
15.1 963 + 3 
24.5 . . 715 + 10 
16.4 1155 +2 
20.2 : 1042 +7 
23.9 1018 +3 
31.6 J . 795 + 15 
30.0 J 685 + 20 
40.5 J 905 + 5 
33.0 J : 1105 + 10 
36.0 ‘ . 1350 + 50 
54.2 720 + ? 


minimum temperature on the dashed curve. No data are avail- 
able on the extent of this two-liquid field into the ternary 
system from the MgO-SiO, side other than that by interpola- 
tion of isotherms from points on the sidelines and composi- 
tions 39, 40, 50, and 80. To these must be added Greig’s 
point at 4.4 per cent K,0, 10.6 per cent MgO, and 85.0 per 
cent SiO,** which had a cristobalite liquidus at 1541 + 9°, 


24 All per cent figures refer to weight per cent. 
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Ww 
J’ 
oO’ 
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V 
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E’ 
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and melted to one liquid (Greig, 1927, p. 40). The above data 
make any extension of the two-liquid field beyond 2 per cent 
K,O improbable, and the agreement of the isotherms drawn 
from the above with the data on the system Na,O-MgO-Si0, 
(Greig, 1927, p. 35) is sufficiently good to warrant the as- 
sumption that the field of two liquids extends about the same 
amount or even less into the K,O system (i.e., to about 1.5 
per cent K,O or less). This is to be expected in view of the 
weakly polarizing nature of the K* ion which aids the forma- 
tion of oxygen bridges and hence tends to eliminate liquid 
immiscibility (Barth and Rosengqvist, 1949). 


Tridymite.—The field of tridymite extends from 1470° down 
to the eutectic point R (NABRQPON, fig. 5) at 715 + 10°. 
Thus it extends metastably over an area where quartz should 
be the stable form. Since K,0 and MgO are not known to 
enter the quartz lattice, but may well be present in appreci- 
able quantities in the tridymite lattice (Lukesh, 1942, found ap- 
proximately 21 per cent CaO, Na,O, and Al,O, in tridymite 
from Plumas Co., this tridymite having the empirical formula 
NaCaAl,Si,;O..), the inversion, under equilibrium conditions, 
should be lowered an appreciable amount below 867° where- 
ever it can occur in the ternary system. No quartz formed 
in any of the mixtures, although a number were crystallized 
at temperatures below the 867° (Kracek, 1939) inversion of 
tridymite to quartz, and two had their metastable tridymite 
liquidus considerably below 867° (110 and 111). Presumably 
these would have a higher temperature stable quartz liquidus 
if held for a longer period, but as these glasses are exceedingly 
viscous they might require excessive time.” 

The difference between the metastable tridymite and the 
higher temperature stable quartz liquidus is unknown in this 
system, but its magnitude may be judged from several ex- 
perimental determinations of it in binary systems. In the sys- 
tem K,0-SiO, the metastable tridymite-K,0°4Si0, eutectic is 
only 2° lower than the stable quartz-K,0°4Si0, eutectic 
(Kracek, Bowen, and Morey, 1937), and in the system Na,O- 
SiO, the metastable eutectic of Na,O°2Si0, with tridymite is 
11° below the stable eutectic with quartz (Kracek, 1930b). 
These small differences in the eutectic temperature may be 


25 Three-month crystallization of 110 at temperatures of 800° and 720° 
failed to yield quartz crystals. 
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misleading, however, as the difference between the stable 
liquidus temperature of quartz and the metastable one of 
tridymite for a given composition is much greater because of 
the steep slope of the silica liquidus and lateral shift in the 
eutectic point. When the data for these two systems are plotted 
on an enlarged scale, the metastable tridymite liquidus is 
found to be lower than that of stable quartz by amounts vary- 
ing from 0° (at the 870° inversion point) to as much as 50° 
for compositions at the metastable eutectic. 

The tridymite crystals formed were exceedingly thin hexa- 
gonal flakes; the size and speed of crystallization varied with 
the temperature of formation. 


Clinoenstatite——The field of the pyroxene clinoenstatite 
(INOTUJI, fig. 5) is unusual in the degree to which it en- 
croaches upon the area where silica might be expected. A com- 
position to the right of N (fig. 5) can have over 82.5 per cent 
SiO, and only 10 per cent MgO in it and still form the pyrox- 
ene clinoenstatite (composition 59.9 per cent SiO, and 40.1 
per cent MgO) as the primary phase. Three of the other 
alkali-alkaline earth-silica systems, Na,O-CaO-SiO, (Morey 
and Bowen, 1925), K,0-CaO-SiO, (Morey, Kracek, and Bowen, 
1980), and Na,O-MgO-SiO, (Botvinkin and Popova, 1937; 
Manuilova, 1937), are similar, but the most siliceous liquid 
from which the pyroxenoid pseudowollastonite (CaO-SiO,) can 
form as the the primary phase is only 75 per cent SiO, in 
the first and 76 per cent in the second; the pyroxene clinoen- 
statite separates out from certain compositions with as much 
as 77 per cent SiO, in the third of these systems. Other ternary 
systems of silica with K,0 and Na,O (Kracek, 1982), with 
diopside and leucite (Schairer and Bowen, 1938), with an- 
orthite and leucite (Schairer and Bowen, 1947b), and with 
forsterite and leucite (Schairer, 1948) show similar tendencies 
but to a lesser degree. 

The field of clinoenstatite is rather narrow, as may be seen 
from figure 5 and the following example: a mixture of the 
composition 68 per cent SiO,, 832 per cent MgO (on the MgO- 
SiO, binary) has a cristobalite liquidus at 1660° (Greig, 
1927, p. 17). The addition of 1 per cent K,O to this composi- 
tion will cause it to form clinoenstatite as the primary phase, 
with a liquidus of 1540°, and the addition of 2.7 per cent more 
K.O (3.7 per cent total) will cause it to form forsterite as 
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the primary phase, with a liquidus of 1533°. Compositions 
near to the clinoenstatite-silica eutectic on the MgO-SiO, 
binary (at 65 per cent SiO,, 35 per cent MgO) need the ad- 
dition of only about 2 per cent K,O to change them from a 
cristobalite to a forsterite liquidus. 

The strong curvature of the clinoenstatite-tridymite bound- 
ary curve (fig. 5) results in a rapid change in the bulk com- 
position of the crystals forming at any given -AT' along this 
boundary curve, e.g., from 87 per cent clinoenstatite, 13 per 
cent tridymite at point N to 13 per cent clinoenstatite, 87 
per cent tridymite just above point O, as shown by the inter- 
section of the tangent to the boundary curve with KS’. 

The crystal habit of the clinoenstatite varied widely. In 
the high-temperature compositions (>1300° +) the crystals 
were short, stubby prisms (almost equant) with the forms 
(100), (010), (110), and (111), an extinction angle (Zac) 
of 22°, the characteristic multiple twinning on (100) (Allen, 
Wright, and Clement, 1906), and pronounced rough cracks 
parallel to (001) (Bowen and Andersen, 1914). At lower 
temperatures (875-1300 +) those compositions which lie 
in the “Clinoenstatite” field of figure 3 crystallized to ex- 
tremely thin needles elongated parallel to the slow vibration 
direction (+ elongation) with parallel extinction, and with 
multiple twinning absent throughout.” In a number of prepa- 
rations these acicular crystals were formed metastably in the 
first few minutes at the start of crystallization, only to 
disappear after a few weeks in the furnace. 

It might appear that this system would offer good pos- 
sibilities of clarifying the confusion that exists concerning 
the various forms of the compound MgO-Si0., but up to the 
present it has not done so. The monoclinic form clinoenstatite 
has been obtained from melts with various fluxes in good 
crystals suitable for goniometry (Allen, Wright and Clement, 
1906). Good crystals of the orthorhombic form enstatite have 
not been obtained in the laboratory.” Grigor’ev (1935) found 

26In the system MgO-Al,O,-SiO, the clinoenstatite crystals formed 


were “seldom twinned” although 1345° was the lowest liquidus temperature 
involved. See Rankin and Merwin, 1918, footnote 2, p. 302. 

27 The second orthorhombic form of MgO-SiO, obtained by Allen e¢ al., 
p. 400 and 433, similar to Kupfferite, was later shown to be a submicroscopic 
intergrowth of forsterite crystals and a siliceous glass, the mass having 


the composition MgO -SiO, and the optical properties of a homogeneous sub- 
stance. See Bowen and Posnjak, 1931. 
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that the addition of CaF, was necessary to cause enstatite to 
form on cooling a mixture of the composition 3MgO:Ca0- 
3Si0,. Bowen and Schairer (1935) caused new crystals of 
orthorhombic (Mg, Fe)O°SiO, to form by the recrystallization 
of natural bronzites, but only with the aid of a small amount 
of NaF. Thilo (1937) converted “metatale” (from the thermal 
decomposition of talc) to enstatite at 790° by the addition of 
LiF. The only syntheses of enstatite not involving the use of 
a flux seem to be formation both from a glass of the com- 
position MgO-SiO, at 1000-1100° (Allen, Wright, and Clem- 
ent, 1906, p. 400) and from thermal decomposition of natural 
anthophyllite and tale (Thilo and Rogge, 1939; Rigterink, 
1947). The former “enstatite” occurred as acicular crystals 
in the glass and had definitely parallel extinction in all orienta- 
tions,” this being the only criterion then available to identify 
it as enstatite. Some of the preparations in the present study 
formed similar length-slow (i.¢., with Z vibration direction 
parallel to the long axis) acicular crystals with definitely paral- 
lel extinction, but as the temperatures were considerably above 
the 1140° inversion, these crystals cannot be enstatite. In 
those cases where thicker crystals grew (mostly above 1800°) 
they usually were of inclined extinction and hence presumably 
clinoenstatite. Bowen and Schairer (1935) give some evidence 
for an enantiotropic relation between enstatite and clinoen- 
statite with an inversion at about 1140°, enstatite being the 
stable form at <1140°. The reaction rate at temperatures 
near the inversion temperature is small. Four hours at 1150° 
were required to effect even a partial transformation. This 
inversion temperature is for the pure compound ; the possibility 
is not excluded that small amounts of K,O may enter the 
pyroxene lattice in the present system and affect the inversion 
temperature slightly. 

Enstatite would not be expected to appear at the liquidus 
surface in the present system on the basis of this inversion to 
clinoenstatite, as all liquids from which MgO-SiO, can crystal- 
lize complete their crystallization above that temperature 
(i.e., at O, 1165 + 2°, or at U, 1155 + 2°, see fig. 5). Under 
equilibrium conditions this inversion should take place in the 
solid state when these compositions are held at temperatures 

28 Allen et al. (1906), have shown that this precaution is necessary. 


due to the occasional tabular habit of clinoenstatite crystals parallel to 
(100) yielding parallel extinction. 
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below 1140°, but no unequivocal evidence was found here to 
prove it, probably due to the slow rate. 

Above 1400° there seems to be little doubt that clinoenstatite 
is the only stable phase, and all other forms, regardless of 
their source, invert to it, the differences being due to different 
starting materials, the rate factor (Allen, Wright, and Clem- 
ent, 1906, above 1300° ; Haraldsen, 1980, above 1400°?; Ewell, 
Bunting, and Geller, 1935, between 1200° and 1300° ; Biissem 
and Schusterius, 1988, above 1400° ; Thilo, 1989, above 1270° ; 
Thilo and Rogge, 1939, above 1270°; and Eitel, 1943, above 
1200°), and impurities (Rigterink, 1947; Haraldsen, 1930, 
reported rather large amounts of impurities in the SiO, and 
MgO used). Below these temperatures there appear to be 
forms of MgO-SiO, other than enstatite and clinoenstatite. 
Haraldsen (1930) found, on the basis of X-ray data, that by 
heating talc a new form of MgO-SiO, developed. Ewell, Bunt- 
ing, and Geller (1935) did not duplicate the results of Harald- 
sen, but Thilo (1939) and Biissem and Schusterius (1938) 
did verify them, and in later work (1939) on anthophyllite and © 
talc Thilo and Rogge discovered two other new forms, “M,” 
and “M,,” in addition to the form found by Haraldsen.” There 
seems to be little doubt that at least some of these forms are 
metastable. Mesoenstatite is not stable at room tempera- 
tures, as Biissem et al. (1939) and later workers (Thurn- 
auer and Rodriguez, 1942) found that grinding to prepare 
samples for X-ray work will cause mesoenstatite (protoen- 
statite) to invert to clinoenstatite unless 25 to 50 per cent of 
glass is present to stabilize it. Clark (1946) reports that 
mesoenstatite inverts spontaneously within two months at room 
temperature to form clinoenstatite with a volume change suf- 
ficient to crumble a fired briquette. Biissem and Schusterius 
(1938) found that various fluxes, such as H,O, Li,O, CaF,, 
H;BO,, and borax, aided the change to clinoenstatite, and 
the inversion seems to be irreversible (Eitel, 1943). 

In the present investigation a number of the samples which 
showed parallel extinction needles were X-rayed by the powder 
technique, but because of the distance between the composition 
points of the mixtures and that of MgO’SiO, the necessarily 
small percentage of crystals present in the glass in each case 


29 The original form found by Haraldsen and the “metatalc” of Thilo 
appear to be identical with the “protoenstatite” of Btissem and of Eitel, 
and the “mesoenstatite” of Thilo and Rogge, and of later workers. 
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yielded rather poor patterns, and unfortunately the patterns 
of all forms of MgO’SiO,—clinoenstatite, enstatite, mesoen- 
statite, M,, and M,—are rather similar in some respects. How- 
ever, several of the patterns obtained were very much more 
like that given by Clark for mesoenstatite than those for 
clinoenstatite or enstatite. A more thorough check on the 
optical properties of the MgO-SiO, crystals obtained in the 
various preparations and more detailed X-ray work are 
planned. Recently, while the present paper was in press, W. R. 
Foster presented a paper on X-ray and thermal data at the 
meeting of the Geological Society of America in Washington, 
D. C., November 1950, indicating that protoenstatite (meso- 
enstatite) is the stable form at moderate and high temperatures 
but that it inverts during the cooling cycle to clinoenstatite 
in a manner similar to the inversion of high to low cristobalite. 


Forsterite——The field of forsterite covers a large part of 
the diagram (LJUYI'G’K’ML, fig. 5). Its extension toward 
the SiO, apex (up to 72.7 per cent SiO, near U, fig. 5) is 
similar to that of clinoenstatite but not as pronounced. The 
isotherms on the forsterite field above 1600° are only ap- 


proximations up to the melting point of forsterite at 1890° 
+ 20° (Bowen and Andersen, 1914), as excessive volatilization 
of the platinum winding at higher temperatures makes 1600° 
the useful working limit of the furnaces. 

The crystals of forsterite formed at high temperatures were 
nearly equant, sharp prisms, usually showing the forms (100), 
(010), (101), (021), (001), and (111). At lower tempera- 
tures (near the line UZI’G’, fig. 5) the forsterite crystals are 
usually very thin plates, blades, or needles, with the Z vibra- 
tion direction parallel to the elongation (length-slow, or 
+ elongation). Although many of these plates and blades were 
elongated to five times their width, the direction of elongation 
was frequently not parallel to any one of the crystallographic 
axes, resulting in inclined extinction of as much as 40°. Al- 
though these crystals were very thin, the elongation appeared 
to be due to excessive development of diagonally opposite 
faces of (011) or (101). 


Periclase——Only one composition was situated on the peri- 
clase liquidus surface (9), although indirect evidence from 
120 and 121 is concordant with the data from it. Both of these 
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latter compositions yielded tiny sharp octahedra of periclase, 
an unusual habit for this mineral since it normally crystallizes 
from melts as completely rounded grains.” X-ray examination 
of this material verified the petrographic examination. In cer- 
tain mixtures in this system periclase will form as the primary 
phase with as much as 45 to 46 per cent SiO, present. Pre- 
sumably there is a smooth rise in the liquidus surface of peri- 
clase up to its melting point at 2800° (Kanolt, 1913; Ruff, 
1913, melting point >2500°). 


Potassium Silicates.—The fields of the potassium silicates 
K,0°4Si0,, K,0°2Si0,, and K,0°SiO, are small, extending out 
into the ternary system to only 214, 6, and 2 per cent MgO 
respectively, because of the low melting points of these three 
silicates and the formation of three ternary compounds with 
a 1:1 ratio of K,0:MgO. Compositions having K,0-4Si0, as 
the primary phase (area CBRD’C’DC, fig. 5) were very re- 
luctant to crystallize, although they eventually formed thin 
blades and needles of K,0°4Si0, characterized by their low 
indices and low birefringence. Compositions having a K,0- 
2Si0, liquidus (area EDC’B’F’M’FE, fig. 5) crystallized 
readily (in a few hours or days) to yield booklets of plates, 
occasionally of parallelogram outline with the optic plane in 
the plane of the tablet and the Z vibration direction bisecting 
the acute angle. Similar crystals have been described by Niggli 
(1913) and by Morey (1914). The K,O0°SiO, field is very nar- 
row and only one composition was made in it. As no composition 
was made with a K,0:SiO, ratio >1, nothing can be said about 
the appearance of the orthosilicate 2K,0°SiO, in the system 
(Kréger and Fingas, 1933). 


Stability Fields of the Ternary Compounds. 

(?). 

The field of this compound occupies a rather large area 
in the lower portion of the diagram (P’K’G’G’F’L’M’'N’, 


80 This rounded growth of crystals seems to be most frequently found 
in very slightly undercooled melts; perhaps the slight oscillation in tem- 
perature due to the normal regulator action is sufficient to result in re- 
peated minor solution and redeposition, with very slight undercooling. The 
preparations yielding the euhedral crystals of periclase were heated in a 
gas-fired furnace, not controlled by a regulator. For a discussion of this 
rounded growth of crystals see Naken (1915, 1917) and Tammann (1925). 
Sharp octahedra of MgO have been obtained from melts containing no 
silica (Rankin and Merwin, 1916). 
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fig. 5). It occurs usually as colorless transparent rounded 
grains, isotropic, with n = approximately 1.54, and a strong, 
distinctive X-ray powder pattern. Many of these grains pos- 
sessed a small gas-filled cavity at the center, and occasionally 
sharp octahedra or skeletal growths consisting of three thin, 
mutually perpendicular dendrites were formed. There were 
some discrepancies in the measurement of the index of refrac- 
tion from sample to sample; these may be due mainly to re- 
action of these high alkali materials with the index liquid, or 
they may be connected with the observed variation in the unit 
cell dimensions (see under X-ray Data). Another possibility 
is that some of these crystals were the compound K,0-MgO 
2Si0, (see footnote 49). 

The evidence is in favor of the composition of most of these 
isotropic crystals being K.0°'MgO’SiO., but is not adequate to 
state so definitely because of the high melting temperatures 
and the high alkali losses in this area of the diagram. When 
a mixture of the composition K,0-MgO'SiO, was made up 
and held at 1620° for 30 minutes,” the batch was found to 
consist of a mass of round, clear, isotropic crystals with a 
small amount of interstitial glass and undissolved MgO. Re- 
melting at 1650° for 10 minutes resulted in fusion and on 
cooling a dendritic mass of the same isotropic crystals with 
more glass was found. MgO was present also but had recrystal- 
lized to the sharp octahedra mentioned under “Periclase.” 
Further fusion of a small sample of this batch for 5 seconds 
in an oxy-gas flame resulted in visible alkali fumes and the ap- 
pearance of an appreciable quantity of forsterite crystals in 
the melt, along with the MgO and isotropic crystals. This be- 
havior is caused by the loss of alkali changing the composition 
of the batch. In one case one 15 minute fusion in air at 1600° 
resulted in a drop of the K,O content from the weighed 48.4 - 
per cent down to 28.1 per cent (determined by weight loss). 
Other attempts to synthesize this compound gave similar re- 
sults. Batch 120, calculated to give approximately 95 per 
cent K,0-MgO:SiO,, 5 per cent MgO, showed the same be- 
havior, but the amount of periclase in the fused material was 
definitely greater than in 121. 

The above data indicate K,0-MgO'SiO, to be the probable 


81 Using a gas-fired furnace and a Leeds and Northrup disappearing 
filament-type optical pyrometer. 
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composition of the compound and 1650° its approximate 
melting point, although the possibility of incongruent melt- 
ing to MgO and liquid is not excluded. More accurate work 
will presumably require the control of the raiher high alkali 
vapor pressures. Several other batches of integral molecular 
composition in the general area (with the ratios 2:1:2, 2:1:3, 
2:1:4, 1:1:2, and 3:2:6 potash to magnesia to silica) were 
made up but cannot be the isotropic compound since they gave 
much lower liquidus determinations with the isotropic crystals 
as the primary phase, wide melting ranges, and definitely 
heterogeneous crystalline products. One of these ratios (2:1:3) 
occurs as an isometric compound in the systems K,O0-CaO-Si0, 
(Morey, Kracek, and Bowen, 1930), Na,O-CaO-SiO, (Morey 
and Bowen, 1925), and possibly K,0-ZnO-SiO, (Ingerson, 
Morey, and Tuttle, 1948), but the data on a batch of that 
composition (119) and neighboring points would seem to pre- 
clude both the possibility that that is the composition of the 
isotropic crystals in the present system and the possibility 
that such a compound occurs above the solidus in this system. 

Apparently the counterpart of K,O-MgO’SiO, does occur 
in the potash-zine oxide-silica system in which an isotropic 
compound tentatively identfied as K,O°ZnO°SiO, was found 
(Ingerson, Morey, and Tuttle, 1948, p. 33). The similarity 
of ionic size of Zn** and Mg*+* (0.88 and 0.78 A) should 
result in similar structures for these two compounds. The sub- 
stitution of the larger Ca** ion (1.06 A) for Mg** cannot 
be made without changing the lattice, however, as the com- 
pound K,0°CaO’SiO, is hexagonal (Morey, Kracek, and 
Bowen, 1930), and the substitution of the still larger Pb* * 
ion (1.82 A) apparently cannot be made at all (Geller and 
Bunting, 1936; McMurdie, 1941). The smaller sodium ion 
(0.98 A) can presumably replace the larger potassium (1.33 
A), as the compound Na,O-MgO°Si0, is isometric also (Bot- 
vinkin and Popova, 1937; Manuilova, 1987) and the substi- 
tution of 2Al*** for Mg** and Sit*** results in the com- 
pound K,0:Al,0, which is also isometric (Brownmiller, 1935), 
and apparently isostructural with K,0-MgO°SiO,. 

K,0°MgO:SiO, forms binary systems with K,0°Si0,, K,0- 
2Si0., K,0°-MgO'3Si0., and probably with MgO (see figs. 5 
and 6). The system with forsterite is only partially binary, 
as the field of periclase covers most of it. 
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K,0°Mg0-3SiO,.—The field of this compound is outlined 
on figure 5 by the line I’Z B’F’G’G'I’. The compound melts 
congruently and sharply to a rather thin liquid at 1134° 
+ 1°, which crystallizes readily upon cooling (within several 
minutes) to a granular mass. Very slow cooling of a liquid of 
this composition through the liquidus temperature yields large, 
sharp hexagonal tablets, uniaxial negative, with n, = 1.530 
+ .002, n, = 1.524 + .002 (Na,**°). In batches of other 
compositions and less often in the pure material, it crystallizes 
occasionally as trillings with hexagonal outline or fourlings 
with a rectangular or circular outline, but these twin crystals 
appeared to be almost optically continuous across the “twin- 
ning planes.” In some cases a series of grooves were seen on 
the (0001) face of the hexagonal tablets, perpendicular to each 
edge exactly as in the pseudo-hexagonal cyclic twins of ara- 
gonite, but such crystals gave sharp, definitely uniaxial inter- 
ference figures and uniform extinction throughout. Several 
runs yielded long, rectangular cross-section crystals of the 
compound with a twinning plane running the length of the 
crystal, the two halves showing 27° extinction against this 
plane. All of the above types of crystals are presumed to be 
the same phase and are here called «K,0-Mg0-3Si0,. This 
is the form which occurred in those quenching runs where suf- 
ficiently large crystals were obtained to differentiate, but 
usually this was not possible. The compound has a definite 
inversion to a finely fibrous appearing, spherulitic form having 
an “extinction cross” parallel to the vibration directions of 
the nichols, usually with the fast ray parallel to the radius of 
the spherulite, although the reverse was found. The gross index 
of refraction of this material lies between the two indices of 
the other form (i.¢., approximately 1.527), and there appears 
to be a fairly large dispersion, but further data on the optical 
properties could not be obtained from these samples. As this 
form is apparently the stable one, it is called “8K,0'MgO- 
3Si0.” here. The differences between the two forms under the 
microscope are not great, particularly in fine materials, and 
hence the inversion was not found until X-ray patterns were 
obtained. The nature of this inversion is unknown, as no 
quenching data could be found to indicate it, and high sensitivity 
differential thermal analysis run on a preparation consisting 
of the 8 form showed no thermal activity from 20° to 1080°. 
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The hexagonal crystals seem to form first in any crystalliza- 
tion, and later slowly invert to the other form, even at tem- 
peratures just below the melting point. Optically the change 
seems to be gradual, but one series of X-ray patterns showed 
a rather abrupt change between samples held for 2 and 3 days 
at 1100°. Density measurements were not made of the a form, 
but there appears to be a volume change upon inversion. For 
additional data on this compound see “X-ray Data” section. 
The compound K,O0°MgO-3Si0, forms binary systems with 
the compounds K,0°2Si0., K,O0°-MgO- 
5Si0,., and forsterite, as shown in figure 6. Glass of the com- 
position K,0-MgO-3Si0, has a density of 2.40 and an index 
of 1.510, and crystals of the 8 form have a density of 2.56. 
K,O-Mg0O-5SiO,.—The field of this compound is outlined by 
the line A’ZX VPQRC’B’A’ on figure 5. This compound melts 
congruently and sharply, though slowly, to a viscous liquid at 
1089° + 1°, and will start to crystallize from a melt of its 
own composition within a few hours at optimum temperatures. 
It occurs in two modifications; an isotropic form called 
“a K,0°-MgO'5Si0,,” as tetrahexahedra and occasional octa- 
hedra and cubes with an index of 1.501 + .002 (Na,*°); 
and a very weakly birefringent (estimated birefringence 0.001- 
0.002) spherulitic, almost fibrous form called “8 K,0-MgO- 
5Si0.,” showing wavy, irregular extinction, and an index of 
1.505 + .002 (Na,**°). The isotropic a form always seems to 
crystallize first, and, if sufficient time is allowed in the furnace, 
it slowly inverts to the anisotropic 8 form, regardless of the 
temperature. The isotropic form has not been obtained free 
from both glass and the anisotropic form; whether this is 
due to the necessity of having glass present to prevent a 
“snap” inversion to the anisotropic form on cooling, or whether 
it is merely a result of inversion proceeding almost as fast as 
the crystallization of the isotropic form, is not known. It is 
very apparent, however, that the glass surrounding crystals 
of the isotropic form, and even that around those crystals that 
have inverted to the anisotropic one, is under strain indicating 
compression along the tangent and tension along the radius 
from each grain. This might be due to a shrinkage of the 
crystal upon inversion as would be indicated by the increase 
in index, but this would not explain similar strains around 
each seemingly perfectly isotropic crystal unless there is a 
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partial inversion that is not visible. Possibly these strains are 
attributable merely to large differences between the thermal 
expansions of the crystalline material and the glass. The 
strains in the glass frequently made it impossible to distinguish 
the two forms in the quenching runs, but in general the an- 
isotropic 8 form seemed far more common, the occurrence of 
the isotropic form being irregular. The inversion with time 
can be followed easily by X-ray methods (see beyond), but 
high sensitivity differential thermal analysis of a preparation 
consisting of the 8 form showed no thermal activity from 20° 
to 1080°. 

K,0°MgO'5SiO, forms binary systems with the compounds 
K,0°4Si0,, SiO,, and forsterite, and partially 
binary systems with K,0°2Si0, and K,0°5MgO'12Si0.,, as in- 
dicated in figure 6. Glass of the composition K,O’'MgO:5Si0, 
has a density of 2.38 and an index of 1.498, and crystals of 
the 8 form have a density of 2.389,. 

K,0°5Mg0O:12Si0O,.—This is the most interesting of the 
new ternary compounds, and was the most difficult to prove 
as to composition because of its incongruent melting and very 
small field on the liquidus surface (see fig. 5, area TOPVXUT). 
Although the field is small, it forms as one of the stable phases 
in all compositions in the triangle LW'S’ (figs. 5 and 8), and 
occasionally is formed metastably in compositions outside of 
that triangle. It generally crystallizes within a few days to 
sharp, thin, hexagonal tablets, frequently with the second- 
order prism and first- and second-order pyramids as modifying 
forms. The small crystals obtained here have no visible cleav- 
age. Optically it is uniaxial positive, n, = 1.543, n, = 1.550 
+ .002 (Na,**°), n,n, = 0.008 + .002. It should be noted 
that these optical properties are very similar to those of 
quartz. Crystals of the compound containing small amounts 
of glass and clinoenstatite (see beyond) had a density of 2.58. 

The thermal behavior of this compound is illustrated by 
figure 7. Upon heating the pure crystals, they melt incon- 
gruently at 1174° + 2° to form clinoenstatite (31.8 per cent) 
and liquid T' (68.7 per cent). Above 1174° this clinoenstatite 
goes into solution, but at 1898° the liquid composition touches 
the forsterite-clinoenstatite reaction curve (JU, figs. 5 and 8) 
and thus the remaining clinoenstatite begins to break up to 
form forsterite and more liquid. At 1453° this reaction is 
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completed with the decomposition of the last trace of clino- 
enstatite; with further heating the remaining forsterite crystals 
start to dissolve and the liquid composition leaves the curve 
JU and approaches point S, where the last of the forsterite 
crystals dissolve at the forsterite liquidus temperature of 
1467°. 

As the composition K,0°5MgO°12Si0, is a little unusual, the 
following discussion is offered in its support. The best method 
of proving the composition of a compound in this type of work 
would be to synthesize an optically homogeneous crystalline 
phase of integral molecular ratio of constituents possessing 
the thermal behavior of a compound. From detailed work 
on the reaction curve OT'U between the fields of clinoenstatite 
and K,0°5Mg0'12Si0, (see fig. 8), it was possible to lo- 
cate the maximum point 7 with fair certainty although the 
temperature differences are small. On a theoretical basis the 
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Fig. 7. Equilibrium diagram for the partially binary system MgO. 
SiO,-K,0.5MgO.12SiO,. This represents a cross section through the 
ternary system along the line KST (fig. 5), extended to its intersection 
with the K,O-SiO, side at the molecular ratio 1K,0:7SiO, (not a com- 
pound). The letters used for the points are the same as those on figures 
5 and 8. Heavy lines refer to binary equilibrium, light lines to non-binary 
(ternary) equilibrium. The small circles are experimentally determined 
points, mostly projected short distances (<0.5 per cent) along the iso- 
therms in the ternary system. 
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composition of the (at that time) unknown phase would thus 
lie along the line connecting T' with the composition of clino- 
enstatite (XK). On the same basis the small maximum at V on 
the boundary curve XVP (between K,0°-MgO-5Si0, and the 
unknown phase) would require that the composition of the un- 
known phase lie on the extension of the line WV. If the maxima 
T and V could be located with complete accuracy, this would 
be adequate to establish the composition of the phase, which 
must lie at the intersection of WV produced and 7K in order 
to satisfy both requirements. Experimentally, however, this 
accuracy was not feasible. When the possible errors in loca- 
tion of T' and V were taken into account a small area in the 
vicinity of point S was found to be the most likely. The only 
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Fig. 8. Equilibrium diagram for the system Forsterite-K,0-.MgO-5SiO,- 
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integral molecular ratio that fell within this area that was at 
all rational was the composition K,0-5Mg0°12Si0,. 

The next step would be to produce a homogeneous crystal- 
line product of this molecular ratio. A number of methods 
were tried unsuccessfully in an attempt to obtain such a pro- 
duct. When a glass of that composition (98 and 98a) was 
held at any temperature below the 1174° incongruent melting 
point, the first crystals to form were metastable “clinoen- 
statite” needles. This occurred within a few hours when held at 
1100° and within a few days at temperatures around 875°. 
Within a day (or more at the lower temperatures) large 
hexagonal plates of the new compound started to grow in the 
glass surrounding the clinoenstatite needles and to enclose 
some of them. With further crystallization of these plates the 
clinoenstatite in the surrounding glass gradually decreased 
in amount until finally a mass of crystals of the new phase 
was obtained, surrounded by films of glass and enclosing a 
small precentage of fine needles of clinoenstatite, after about 
three days at 1100°. Actually this glass was of composition 
T (figs. 5 and 8) and was in the stoichiometrically correct 
amount to combine with all of the clinoenstatite remaining to 
form more of the new phase, but as the glass-new phase inter- 
face and the new phase-clinoenstatite interface are each in 
themselves stable at these temperatures, the clinoenstatite is 
effectively removed from further reaction with the glass and 
solid diffusion through the lattice of the new phase is the only 
way that the reaction can continue. The rate of this diffusion 
at 1100° is apparently slow enough to be experimentally 
negligible, as an additional 17 days at 1100° failed to produce 
any further visible change, particles of the new phase only a 
few microns in diameter still possessing their films of glass 
and enclosed clinoenstatite needles. Fluxes such as potassium 
tungstate were tried, but unsuccessfully. Even when a glass 
of the composition K,0°-5Mg0-12Si0, was carefully “seeded” 
before crystallization with crystals of the compound containing 
no clinoenstatite, the metastable clinoenstatite still formed. 
Since the composition of the glass cannot be independently 
determined with any accuracy, and the relative amounts of 
glass and metastable clinoenstatite cannot be measured with 
any certainty in material of this type, the composition of 
the compound has not been proven by these data. 
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In order to prevent this formation of metastable clinoen- 
statite, crystallization of glasses of the composition K,0- 
5MgO'12Si0, (98 and 98a) was attempted at 500° for 14 
hours under a pressure of 6000 psi water in a bomb of the 
Morey type (Morey and Ingerson, 1937). The result was a 
mixture of two phases, containing approximately 10 per cent 
of large crystals of a hydrous prismatic material, and 90 per 
cent of a hydrous, fibrous, highly birefringent material with 
indices approximately 1.565 and 1.585. Both of these de- 
composed to an isotropic, low index material on being heated 
to 950° for three minutes. When this heated material was held 
in air at 1100° for one day, it yielded a mass of metastable 
clinoenstatite needles in glass. 

As higher temperatures tend to prevent the formation of 
hydrous phases in bomb runs, samples were run at 900° and 
15,000 psi water pressure for 1 hour in O. F. Tuttle’s (1948) 
hydrothermal quenching apparatus. These runs were success- 
ful, yielding a product which was a mass of interlocking 
crystals with optical properties and X-ray powder patterns 
identical with those of the crystals formed in the “dry” manner 
but without the clinoenstatite. A small fraction of 1 per cent 
of some low index, extraneous material was present, presum- 
ably because slight impurities in the starting material (par- 
ticularly the MgO) caused divergence of the batch from the 
theoretical composition. If the actual composition of the com- 
pound were such that it differed from 98 in having only 0.5 
per cent more MgO, the crystalline material obtained here 
from 98 should contain only 94.5 per cent crystals of the 
compound, the other 5.5 per cent being either glass or crystals 
of K,0°'MgO'5Si0, and SiO,. This amount of such low index 
materials should be visible readily under the microscope when 
immersed in a liquid having the mean index of K,0°5Mg0: 
12Si0, crystals; if the compound contained 0.5 per cent less 
MgO than 98, the crystalline material obtained should consist 
of approximately 97.5 per cent of the compound and 2.5 per 
cent clinoenstatite, and the latter should be even more visible 
under the microscope as a result of its high index. Furthermore, 
even these differences in composition from that of 98 would 
not be sufficient to yield any rational formulas for the com- 
position of the compound. The closest integral ratio of any 
simplicity and probability according to present knowledge of 
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silicate structure” is K,0-4MgO-10Si0,, but this composition 
is so far removed from 98 (it would have 11.0 per cent K,O, 
18.8 per cent MgO, and 70.2 per cent SiO,) that, although 
the quenching data might possibly be interpreted in this light, 
the crystallization of 98 to a homogeneous crystalline phase 
prohibits it. 

In order to prove that there was not a range of compositions 
that would have behavior similar to that of 98 upon crystal- 
lization, two other batches were given the same treatment at 
900° and 15,000 psi water pressure. They were 17 (approxi- 
mately 3.2 per cent more MgO) and 134 (approximately 3.8 
per cent less MgO). The former yielded a mass of K,0-5MgO- 
12Si0, crystals containing a mat of tiny clinoenstatite needles 
(theoretically 16 per cent by weight), and the latter a mass of 
K,0°5MgO°12Si0, crystals coated with films of a low index 
material, presumably K,O0-MgO-5Si0, and SiO, or glass.** The 
above data are considered adequate to establish with fair cer- 
tainty that K,O°'5Mg0-12Si0, is the composition of the com- 
pound. It is stable in equilibrium with K,0°-MgO°5Si0,, forster- 
ite, clinoenstatite, and silica, although most of these sub-systems 
are not binary at liquidus temperatures due to the incongruent 
melting of the compound. 


82 E.g., see Hendricks (1944). A series of infinite sheets of (Si,O,,),-4" 
as in the micas, bonded with K and Mg, would satisfy either K,0-4MgO. 
10SiO, or K,O.5MgO-12Si0,, but layers of 12 coordination. K*+ cannot 
exist as they do in the micas, as the crushed crystals showed no sign of 
the perfect cleavage of the micas. 

33 X-ray examination of these examples was not very satisfactory due 


to the overlap of strong K,O-5Mg0-12Si0, lines upon weak clinoenstatite 
and K,O.Mg0O.5SiO, lines. 


(To be continued) 
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EXPERIMENTS ON THE MOVEMENT 
OF SHELLS BY WATER’ 


HENRY W. MENARD anp A. J. BOUCOT 


ABSTRACT. Experiments in a laboratory flume indicate that the current 
velocity necessary to move terebratuloid brachiopod shells is only about one 
tenth of the velocity required to move pebbles of the same size. Movement 
of shells is controlled primarily by the relative ease of movement of the 
subjacent bed material; a shell may move over a bed of gravel but sink 


into a bed of sand because of the rapid removal of adjacent sand grains 
by scouring. 


Shape, size, effective density, and ornamentation influence the competent 
velocity of a shell. Competent velocity increases with sphericity. 


INTRODUCTION 


ALEOECOLOGISTS are aware that a fossil assemblage 

of shells does not necessarily represent a community as- 
semblage preserved in place. Instead, the fossil assemblage may 
constitute a diverse assortment of shells brought together by 
physical agencies, ordinarily after death, and then buried. 
Moving water is one of the most important agents affecting 
the transportation and sorting of unattached dead shells or 
those otherwise incapable of independent movement. The 
writers know of no experimental observations on the current 
velocities required to initiate the movement of shells, nor of 
the effect of transportation on an assemblage. For these rea- 
sons it seemed desirable to make such observations when a 
hydraulic flume temporarily became available at the Woods 
Hole Oceanographic Institution in the fall of 1948. 

During the experiments, the relations between the variables 
affecting the movement of the shells seemed to be simple; so no 
exhaustive study was attempted. Subsequent analysis of the 
observations revealed that the relations were unexpectedly 
complex. It is apparent now that many more experiments will 
be required before the initial movement of shells can be ex- 
plained completely or predicted. Consequently, the present con- 
clusions cannot be offered as more than first approximations. 
They are submitted at this time because neither of the writers 
is likely to have access to a hydraulic flume in the near future. 

This paper presents some results of the tests, as well as a 
brief consideration of their significance with regard to paleo- 
ecology. The writers wish to express their appreciation to 


1 Contribution No. 546 from the Woods Hole Oceanographic Institution. 
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P. E. Cloud, W. C. Krumbein, Henry C. Stetson, Fred B. 
Phleger and J. Wyatt Durham for reading the manuscript 
and offering suggestions for its improvement. Dr. Cloud initi- 
ated the statistical analysis of the data. In this analysis he 
was advised by F. M. Wadley; W. Peddle made most of these 
computations. The writers are also indebted to the Woods 
Hole Oceanographic Institution for the use of laboratory 
facilities, particularly the flume. 

The experiments were done in a wooden flume one foot wide 
and 6 feet long, with a glass observation panel in one side. In 
26 tests, a recess in the central 3 feet of the flume was filled 
with medium-fine sand (diameter 0.36 mm.) to the level of the 
flume bottom; in one test, the recess was filled with a gravel 
composed of grains 2-5 mm. in diameter. The bottom of the 
flume had a fixed slope of 0.002 during all these tests. In a 
typical run, a tailgate was fixed in place in order to produce 
the desired depth of flow in the flume; several shells were 
placed on the smoothed sand bottom; and the flume was filled 
with fresh water at a rate slow enough so that neither the 
shells nor the sand were disturbed; then the flow of water in 
the flume gradually was increased until one of the shells or 
the sand moved. The depth of water in the flume and the dis- 
charge of the water were measured when movement occurred. 
The mean current velocity was then calculated by dividing the 
discharge by the cross-sectional area of the flow within 
the flume. 

Some 27 runs were executed, and 39 individual dry shells, 
representing 9 species of terebratuloid brachiopods, were 
tested. Observations were made on the current velocity com- 
petent to move the shells, the nature of the movement, scour 
of sand around the shells, and the orientation of the shells dur- 
ing and after movement. The length, width, thickness, weight, 
nominal diameter, and sphericity of the shells used in these runs, 
and the results of the tests are shown in tables 1 and 2. 

The velocity required to move a brachiopod shell lying free 
on a substratum of sand depends on the shell properties of 
size, effective density, shape, and ornamentation, and on 
whether the valves are hinged together or separated. All these 
variables except ornamentation are considered in the present 
tests. The most important variable governing the initial move- 
ment of shells, however, appears to be the relative ease of 
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movement of the shell and the underlying sand grains. Thus, a 
certain current may roll a shell over a gravel bottom with ease, 
but be incapable of moving it over sand because the sand 
moves more readily than the shell. If the sand moves first, a 
scour is formed around the shell, which may be dropped into 
it and buried by moving sand. Because of the importance of the 
relative ease of motion of sand and shells, the detailed discus- 
sion of the tests is divided into two sections, one concerned 
with shells which move more easily than sand, and the other 
concerned with shells which move less easily than sand. 


WHOLE SHELLS WHICH ARE EASIER TO MOVE THAN 
THE UNDERLYING SAND 


Unattached shells are clastic particles in the same sense as 
grains of sand or gravel, and hence the factors governing their 
movement are the same. As the movement of sand and gravel 
has been studied in great detail whereas the movement of shells 
has not, it seems profitable to consider sediment movement 
before turning to shell movement. 

The characteristics of a grain which influence its initial 
movement on the bed of a watercourse are the size, density, and 
shape. Hydraulic engineers generally agree that the current 
velocity competent to move coarse sand and pebbles is pro- 
portional to the square root of the diameter of the grains 
(Hjulstrom, 1935); the same proportionality probably holds 
for all larger sizes of gravel (Menard, 1949). The competent 
velocity commonly is thought to be proportional to the square 
root of the effective density (weight per unit volume if the 
grains are immersed in water). As to the shape, Krumbein 
(1942) finds that the more spherical a grain of gravel, the 
slower the current required to put it in motion. 

The present tests show that the initial movement of tere- 
bratuloid brachiopods depends primarily on the effective 
density because that variable appears to govern whether a shell 
moves more easily than the sand on which it rests. Secondarily, 
the velocity required to initiate motion appears to depend on 
the shape and size of the shell. Although the movement of 
shells and of gravel and sand is necessarily controlled by the 
same physical characteristics, it is apparent that the relative 
importance of these characteristics varies. 
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EFFECTIVE DENSITY 


All the whole Terebratulina shells tested were moved by cur- 
rents which were incapable of moving the bed of sand, although 
the shells weighed as much as 0.605 gram and the sand grains 
weighed only a few thousandths of a gram. This seeming para- 
dox is due to the very small density of the shell and its con- 
tained water as compared to the solid quartz grains. If shell 
number 4, for example, is filled with water, the shell-water com- 
bination has a density in air of 1.28 gms./cc.; however, if the 
combination is under water, the buoyancy of the displaced 
water reduces the effective density to 0.28 gm./cc. A water- 
filled shell moved by a slow current is analogous to an inflated 
balloon moved by the faintest breeze although the rubber in 
the balloon weighs several grams. 

Species and individuals differ in effective density, and, there- 
fore, in competent velocity. Hence moving water may tend to 
sort a mixture of individuals of two species into groups ac- 
cording to species—rather than size or shape—if the effective 
densities of the two species are very different and the range 
in size and shape is small. Some effective density determina- 
tions for terebratuloid brachiopods appear in table 1. The 
values show that the effective density of Terebratulina septen- 
trionalis ranges from 0.28 to 0.43 gm./cc., and that it de- 
creases as the size increases. The effective density of Tere- 
bratalia transversa ranges from 0.35 to 0.42 but does not vary 
systematically with size. 

Studies of sediment movement indicate that competent ve- 
locity is proportional to the square root of effective density. 
Applying this proportion to brachiopod shells suggests that 
the competent velocity of a one-centimeter 7’. transversa shell 
is only slightly larger than that of a T'. septentrionalis shell of 
about the same size and shape. 


The size of irregularly shaped shells may be expressed in 


terms of a nominal diameter equal to that of a sphere of the 


1A dry shell was weighed in air, and the approximate volume was 
determined by dividing the weight by the density of calcite (2.72 gms./cc.). 
The shell was then filled with water and weighed. The weight of the dry 
shell was substracted from the weight of the water-filled shell in order 
to obtain the weight of the water, and thereby its volume. The density of 
the water-filled shell was calculated by dividing the sum of the volume of 
the shell and of the contained water by the total weight. 
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same volume. For the present purposes, a nominal diameter 
has been determined by taking the cube root of the product of 
the length, width, and thickness of each shell. Thirty-one 
shells (including single pedicle and brachial valves as well as 
whole shells) with nominal diameters ranging from 0.41 to 
5.50 millimeters were moved while resting on a substratum of 
either medium-fine sand or gravel. The current velocity re- 
quired to initiate the motion of each of the shells is shown in 
graphic form in figure 1, in which the competent velocity in 
centimeters per second is plotted against the nominal diameter 
in centimeters. The graph paper is logarithmic because ex- 
periments with sediments have shown that functions of these 
two variables yield a straight line on logarithmic paper. This 
relation provides an easy means of establishing whether or not 
the initial movement of shells is controlled by the size of the 
shells in the same way that the initial movement of sediment 
is controlled by the size of the grains. 

The points relating competent velocity and nominal dia- 
meter in figure 1 are scattered, but the area in which the points 
lie has a trend which suggests that the velocity is proportional 
to the square root of the diameter, as it is for sand and gravel. 
The coefficient of correlation for a least squares regression 
line based on all the data in which the bed was sand gives a low 
positive correlation of 0.585.? If the data for a gravel bed 
are included, a high positive correlation of 0.820? is indicated. 
However, combining data from different experimental condi- 
tions introduces some bias into the calculations. The scatter- 
ing of points indicates that the nominal diameter of a shell 
is not the only variable affecting its competent velocity. This 
indication agrees with the general body of hydrodynamical 
theory and observation. Rouse (1938, p. 216) says that if 
the wake behind an object is completely turbulent, the drag 
exerted upon the object by a moving fluid is almost completely 
independent of the Reynolds number (a criterion of turbulence). 
The Reynolds number in the present experiments is the product 
of the current velocity and the nominal diameter divided by 
the kinematic viscosity; the value of the Reynolds number is 
great enough to establish that the wake was always turbulent. 
With regard to shells (nominal diameter larger than about 
one centimeter), therefore, Rouse’s statement indicates that 


2Computed by W. Peddle. 
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the current velocity necessary to initiate motion may even be 
independent of the size. 


SHAPE 


The shape of shells may be expressed in terms of true 
sphericity—the surface area of a sphere divided by the sur- 
face area of a shell of the same volume (Wadell, 1932). Be- 
cause the surface area of a shell is difficult to measure, the 
shape is here expressed in terms of ratios of diameters—the 
intercept sphericity of Krumbein (1941). Intercept sphericity 
is equal to the cube root of the following quantity: the inter- 
mediate diameter times the smallest diameter divided by the 
square of the largest diameter. This is an approximation of 
sphericity as defined by Wadell. In the present tests, the thick- 
ness is always the small diameter, but either the width or 
length may be the intermediate or large diameter. Sphericity 
for the 31 shells which were moved had a range of 0.522 
to 0.859. 

It will be shown in the following paragraphs that there is a 
correlation between the intercept sphericity of a shell and the 

NOMINAL DIAMETER CM. 


COMPETENT VELOCITY CmM/SEC. 


i 

4 8 
SPHERICITY 
Fig. 1. Competent velocity of terebratuloid brachiopod shells in relation 


to nominal diameter and sphericity. The competent velocity of pebbles is 
shown for comparison. 
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current velocity required to put it in motion. The correlation 
appears to indicate that intercept sphericity is an index of the 
shape of terebratuloid brachiopod shells; but the correlation 
does not mean that the intercept sphericity represents a com- 
plete expression of the shape of the shells. Intercept sphericity 
represents an approximation of true sphericity which holds be- 
cause of the general resemblance between pebbles and a triaxial 
ellipsoid. If the intercept sphericity is a valid expression for 
the shape of only those shells which resemble triaxial ellipsoids, 
it will prove of little value in measuririg the shape of most 
shells which, unlike terebratuloid brachiopods, are irregular. 
It is interesting to note, however, that in the present experi- 
ments with single valves the intercept sphericity appears to 
give a more accurate expression of shell shape than does the 
true sphericity. True sphericity appears to have little signifi- 
cance as a measure of the shape of single valves because each 
valve is essentially a curved plate, and this sphericity takes no 
account of the curvature of such a plate. 

Engineers divide the total drag exerted on an object by a 
moving fluid into two parts: “surface drag” resulting from 
the retardation of the fluid at the solid boundary, and “form 
drag” resulting in large part from the pull exerted by the low 
pressure area in the wake of the particle. Neither true 
sphericity nor intercept sphericity is an exact expression for 
the sum of surface and form drag; hence neither is capable of 
measuring the total effect of shape upon the hydraulic be- 
havior of a particle. True sphericity appears to be an index 
of surface drag but it takes little account of form drag. Con- 
sequently it gives an indication of the hydraulic behavior of 
particles of approximately the same shape, but not of particles 
with very different shapes, because the wakes are different. 
Intercept sphericity gives some measure of the surface drag 
of ellipsoidal particles and also of form drag, because it is 
based on ratios of diameters. With regard to hydraulic be- 
havior, therefore, intercept sphericity may surpass true sphe- 
ricity as an index of shape of sand and gravel particles. 
Irregular shapes such as shells are far more complex, and their 
hydraulic behavior will not be understood completely until the 
principles of fluid dynamics are applied to their study. Neither 
true nor intercept sphericity, for example, is capable of in- 
dicating any significant difference in the hydraulic behavior 
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of a rectangle and of an airfoil of about the same dimensions, 
although the form drag of the airfoil is much less than that 
of the rectangle, and the lift resulting from the shape of the 
airfoil might have a pronounced effect on its movement. A sec- 
tion of a pecten perpendicular to the hinge line resembles an 
airfoil, so it is doubtful if any type of sphericity would provide 
a significant index of its shape in relation to hydraulic behavior. 

Figure 2 shows the sphericity of the various shells graphed 
against the current velocity required to put them in motion. 
For purposes of comparison, the figure also shows the effect 
of particle sphericity on the settling velocity and competent 
velocity of synthetic pebbles with a nominal diameter of about 
14.5 mm. (Krumbein, 1942), and also on the settling velocity 
of molded lead particles (Wadell, 19384). The points relating 
sphericity and competent velocity of shells are scattered, and 
the trend of the relation is subject to several interpretations. 
The four most reasonable interpretations are shown in figure 
3. In the first and second interpretations, the data for single 
valves and whole shells are combined; therefore the relation 
between sphericity and competent velocity is assumed to be con- 
tinuous as shown by Wadell (1934). In the third and fourth 
interpretations, the data for single valves and whole shells 
are separated ; therefore it is assumed that the relation is dis- 
continuous (but continuous for various classes of shapes, such 
as rollers, disks, or spheriods) as shown by Krumbein (1942). 
In the first and third interpretations, the data for movement 
on a sand bed are combined with the data for movement on a 
gravel bed; in the second and fourth interpretations, these 
data are not combined. Statistically there is reason to separate 
the data concerning movement on two sizes of bed material 
because the experimental conditions were different ; combining 
the data may bias any interpretation. Hydraulically the bias 
does not appear to be important. The change in bed type 
probably altered the boundary conditions of the flow in the 
flume; this probably altered the distribution of velocity in the 
main body of the flow but the average alteration probably was 
minor. The shells were all large compared to the thickness of 
any laminar boundary layer, and therefore most of the surface 
of shells on either the sand or gravel bed was exposed to the 
relatively undisturbed turbulent flow outside the boundary 
layer. The data are subject to too many interpretations to 
justify any more extensive discussion of this bias. 
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Despite their diversity, the interpretations indicate that 
there is some correlation between the sphericity of a shell and 
the current velocity required to move it. Some measure of the 
correlation can be determined by drawing a least square re- 
gression line for each interpretation, and calculating a co- 
efficient of correlation for the deviation of the points from that 
line. Interpretation one (figure 3) gives a positive correlation of 
0.699* for the straight line; no coefficient of correlation can be 
calculated for a curved line. “Positive” applied to a correlation 
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Fig. 2. Effect of sphericity upon the competent velocity and settling 
velocity of various kinds of particles. 
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coefficient refers only to its algebraic sign. It is not an indica- 
tion of certainty. Interpretation two gives a low positive cor- 
relation of 0.603 for whole shells.* Interpretation three gives 
a low positive correlation of 0.588. Interpretation four gives 
a positive correlation coefficient of 0.069 for whole shells, in- 
dicating no correlation; for all single valves together it gives 
a low positive correlation of 0.495; for single valves of one 
species, Terebratulina septentrionalis, it gives a positive cor- 
relation of 0.659. Thus almost all of the interpretations suggest 
that the competent velocity of a terebratuloid brachiopod 
shell depends on the sphericity of the shell; but part of the 
fourth interpretation indicates no such dependence for the 
limited number (16) of observations of the movement of whole 
shells included in the interpretation. 

Sphericity is related to shape insofar as the hydraulic be- 
havior of these shells is concerned, but the relation is only ap- 
proximate. Hence two shells with the same weight, size and 
sphericity may have different competent velocities because of 
the shape effect which is not included in the sphericity. This un- 
measured shape effect may account for some of the scattering 
of points in the graphs relating sphericity and competent 
velocity. Scattering may also be due to the action of turbulent 
eddies which cause the instantaneous velocity around a shell! to 
vary from the mean velocity. 

The present tests with terebratuloid brachiopod shells rest- 
ing on a substratum of sand indicate that the more spherical 
shells require a faster current to put them in motion. This 
qualitative relation between sphericity and competent velocity 
is diametrically opposed to a relation observed by Krumbein 
(1942). Krumbein molded a number of artificial pebbles which 
had about the same weight and nominal diameter but different 
sphericities. He placed these pebbles in a flume and then ob- 
served the velocity required to initiate motion as well as the 
relative velocity of the current and the pebbles once they began 
to move. The competent velocities were not published ; however, 
in a letter to Menard in 1949, Dr. Krumbein was kind enough 
to state that the competent velocity for spheres averaged 
14.5 cm./sec., and that for pebbles with a sphericity of about 
0.5 the competent velocity averaged 42.2 cm./sec. Thus the 
velocity required to initiate motion appeared to vary inversely 

+Computed by W. Peddle. 
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with the sphericity. Computations by the present writers using 
Krumbein’s published data on critical Froude numbers yield 
the competence-sphericity relation shown in figure 2. 

On the other hand, Krumbein (1942) also measured the set- 
tling velocity of the artificial pebbles, and his figure 4 (p. 627) 
shows that the settling velocity varies directly with the sphe- 
ricity. Wadell (1934) had noted a similar relation between 
settling velocity and sphericity. Both Krumbein’s and Wadell’s 
data are plotted in figure 2. The settling velocity of particles 
larger than a few millimeters in diameter is directly propor- 
tional to the competent velocity (Hjulstrom, 1935; Menard, 
1949). Consequently, with regard to the relation between com- 
petent velocity and sphericity, Krumbein’s tests with settling 
velocities agree with the flume tests of the present writers. 

The contradictions in observations of the competence- 
sphericity relation of particles in flumes may be more apparent 
than real. A particle is put in motion in a flume if the propelling 
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force due to moving water exceeds the retarding force due to 
friction between the particle and the adjacent particles of 
the bed. Particle shape may have some influence on the effective- 
ness of both the propelling and the retarding forces. Tenta- 
tively, the present writers believe that the differences in the 
observed influence of sphericity on competence may be ascribed 
to the differences in the effective density of the particles in the 
two experiments. The effective density of the shells is less 
than one-fifth of the effective density of Krumbein’s synthetic 
pebbles. The shape of the shells may have more influence on 
the effectiveness of the drag due to the moving water than on 
the effectiveness of the friction between the shells and the sand 
bed. As a result the relation between sphericity and competence 
of shells in a flume corresponds to the relation between sphe- 
ricity and settling velocity (in which there is no friction with 
adjacent particles). With denser particles, however, the pro- 
portion of the grain surface exposed to moving water is much 
smaller compared to the surface in contact with other particles. 
Therefore, the influence of particle shape on the effectiveness 
of the retarding force due to friction between grains may be 
more pronounced than the influence on the effectiveness of the 
propelling force due to’ moving water. If so, the relation be- 
tween sphericity and competence of solid particles in flumes 
would depend on the relation between sphericity and friction 
between grains, and could be just the opposite of the relation 
between sphericity and the settling velocity of the particles. 
Granting the validity of both sets of apparently contradictory 
observations, the fluid velocity required to move shells probably 
varies directly with the sphericity, but the fluid velocity re- 
quired to move most solid detritus (density about 2.7 grams 
per ce.) probably varies inversely with the sphericity. 
Although the conclusion is speculative, some field observa- 
tions also suggest that the velocity required to move shells 
varies directly with the sphericity. DuBois (1916) found that 
individuals of Terebratalia obsoleta (Sowerby) taken from 
Puget Sound varied in shape according to the degree of move- 
ment of the waters in which they lived. The more spherical 
shells were developed in the fastest currents. As the variation 
in shape represents an adaptation to differences in current 
velocity, the more spherical forms must be those which are 
best adapted to live where the currents are faster. The shape 
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may be related to some secondary characteristic of fast cur- 
rents—such as an increased oxygen or food supply—but if 
the adaptation of shape is necessary in order to counterbalance 
the velocity of the water, a more spherical shape must be the 
one which is least affected by moving water. Therefore, the 


velocity necessary to initiate motion appears to vary directly 
with the sphericity. 


ORNAMENTATION 


The effect of variations in ornamentation upon the movement 
of shells can not be evaluated from the present tests because 
all the shells are relatively smooth. Variations in ornamenta- 
tion, however, are variations in shape, so it is apparent that 


they influence the relative ease with which a shell may be put 
in motion. 


HYDRAULIC FACTORS 


The initial movement of a shell on the sandy bed of a 
watercourse is influenced not only by the characteristics of 
the sand and shell but also by the flow characteristics of the 
running water. Friction retards the layer of water next to 
the bottom and sides of the channel so that the velocity of the 
water on a sand bed may be only a fraction of the mean 
velocity of the current. However, the retarded or “boundary” 
layer rarely is more than a millimeter thick ; all but the smallest 
shells protrude through it and are exposed to the force of the 
current. Medium and fine sand grains, however, usually lie 
within the boundary layer and thus may remain undisturbed 
on a stream bed although shells which are larger and heavier 
are in motion. 

The measured effective densities of terebratuloid brachiopods 
are so small that the initial motion of any particular shell 
probably is influenced by the degree of turbulence of the mov- 
ing water. Hjulstrom (1935) has shown that all natural 
watercourses are characterized by turbulent, eddying motion. 
The eddies may move in any direction relative to the general 
motion of the water, with the result that at any instant some 
parcels of water move faster than the average and some move 
slower. Due to the entirely random motion in a turbulent fluid, 
a comparatively rapidly moving eddy may strike a large shell 
and put it in motion although a smaller shell only a few centi- 
meters away remains fixed in position because the force of the 
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water in that place is much smaller. The analogous effect of 
turbulence in air may be observed when a gust of wind swirls 
a patch of dust upward while the surrounding dust remains 
in place, or when a “cats-paw” ripples a water surface in 
some spots but leaves the remainder smooth. These common- 
place phenomena indicate how turbulence may cause the larger 
of two shells of the same species to be moved instead of the 
smaller. 

The quantitative importance of the turbulent eddies with 
regard to initiating motion in whole shells depends on the 
comparative velocity of the eddies and the current as a whole. 
This relation has been the subject of a great deal of study by 
hydraulic engineers, and it is knova that the importance of the 
eddies increases with the cross-sectional area of the current if 
all other things are equal. In a small stream, for example, the 
eddy velocity might be only a fraction of the current velocity, 
but in a large slow current on a continental shelf the eddies 
might flow several times as fast as the average current velocity. 

Perhaps the most significant conclusion to be derived from 
this analysis is that some terebratuloid brachiopod shells may 
be moved so readily that they are capable of being transported 
by the low current velocities now found on the continental 
shelves and slopes. All sizes of Terebratulina septentrionalis 
tested were moved by currents with velocities of 11.7 to 15.1 
em./sec.; Shepard (1948) has recorded bottom velocities of 
36.7 cm./sec. in 160 fathoms on the California continental 
borderland, and 14.6 cm./sec. in 1050 fathoms in the Santa 
Cruz basin. These velocities were measured in the open ocean 
where tidal effects were negligible. However, in bays and 
estuaries where the tidal range is large, current velocities 
commonly are several times as great. 


WHOLE SHELLS WHICH ARE MORE DIFFICULT TO MOVE THAN 
THE UNDERLYING SAND 


Flume tests with individuals of Terebratalia transversa 
show that the transportation of shells may depend on a 
critical relation between the competent velocity of the shells 
and of the sand in the bed. The smallest of five individuals 
resting on a bed of monodisperse sand grains with a median 
diameter of 0.36 mm., was moved by a velocity of 15.0 cm./sec. 
At the same time the sand was scoured from around the re- 
maining four larger shells. A faster current caused the whole 


| 


Experiments on the Movement of Shells by Water 147 


surface of the bed to move although the shells remained fixed. 
They then began to settle into holes in the sand which were 
produced by the scour around them. This sequence of events 
is depicted in figure 4. Similar scour around pebbles was noted 
in earlier experiments in the flume, and fossil scour marks have 
been observed by Peabody (1947) in the Moencopi formation of 
the Grand Canyon area. If a current moves the subjacent sand 
instead of a shell the shell may become buried. However, this 
burial may depend on whether the sand is being eroded or 


Fig. 4. A terebratuloid brachiopod (Terebratulina septentrionalis) being 
buried in sand because of scour by moving water. Four stages in the burial, 
of which the first is uppermost, are depicted in cross sections parallel to 
the direction of flow. In all the stages the water is moving from right to 
left, and at a constant velocity. 
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deposited; if erosion prevails, the shell remains exposed, 
whereas if deposition prevails, it is covered. 

It follows that the size of the grains in the bed may deter- 
mine whether a shell is moved by water with a given current 
velocity. A current which could move a shell on a bed of coarse 
sand might cause it to be buried in a finer sand, because the 
latter would be scoured from around the shell. 

The effect of the size of particles in a substratum is shown 
by the tests with shells number 31 and 39. If these shells were 
resting on medium sand, they began to be buried when the cur- 
rent velocity was less than 20 cm./sec. and they did not move 
even when the velocity was increased to the maximum that the 
flume could produce—about 35 cm./sec. On the other hand, 
both shells were put in motion by a current with a velocity of 
26.8 cm./sec. if they were resting on a substratum of gravel. 

In order to establish whether an assemblage of fossil shells 
can have been transported over a bed of the sand in which it 
is enclosed, it is necessary to know the competent velocity re- 
quired to move both the sand and the shells. The competent 
velocity for any size of sand is known within certain limits 
(Hjulstrom, 1935), but the velocity required to move any 
particular shell can be found only by the type of experiment 
reported in this paper. 

SINGLE VALVES 


Five pedicle valves of Terebratulina septentrionalis, weigh- 
ing 0.055 to 0.407 gm. and oriented concave upward, were 
tested simultaneously. All were moved small distances by a 
current with a velocity of only 11.1 cm./sec. In other tests 
these shells were oriented convex upward, in which position 
one valve weighing 0.115 gm. was moved by a current of 11.1 
cm./sec., and the other valves were moved small distances by 
a current with a velocity of 14.9 cm./sec. 

Another series of tests featured the six brachial valves of 
the same species. These valves weighed 0.042 to 0.257 gms. 
Their behavior at several current velocities is shown in table 2. 
If the current velocity was increased to more than 14.4 cm./sec. 
the sand bed began to move and ripples were formed. The valves 
(convex upward) continued to be fixed in position until a sand 
ripple moved past. Each then was turned over and moved into 
the trough behind the ripple; it remained in the trough as the 
ripple slowly moved downstream. 
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MOVEMENT OF SHELLS 


Most of the whole hinged shells moved over the sand bed by 
rolling with their longest dimension perpendicular to the 
current; however, others slid or rolled end over end. At higher 
current velocities almost all of these latter shells would prob- 
ably roll with their longest dimension perpendicular to the 
current, because this position affords the least resistance to 
movement. Most of the sliding shells maintained a constant 
orientation relative to the current, but the individuals usually 
had different orientations. Both sliding and rolling were epi- 
sodic; a shell slid a few centimeters in a few seconds, then 
remained fixed for several minutes, and slid or rolled again. 
The episodic nature of the movement may have been due to 
the intermittent action of eddies which locally increased the 
mean velocity enough to initiate movement, but soon dissipated 
and no longer moved the shell. 


The Terebratulina shells moving over an immobile sand bed 
came to rest at random orientations relative to the direction 
in which the current was moving. In a fossil state these shells 
would give no clue to the direction of the water movement at 


TasLe 2 


Movement of Brachial Valves of Terebratulina septentrionalis 
(Initial Orientation Concave Upward) 


ton) Behavior At Current Velocity (cm./sec.) 
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the time they were deposited. What is more significant to 
paleoecology is that a small number of shells may have a 
chance parallel orientation, but that there may be any relation 
between the shell orientation end the current direction. 
Single brachiopod valves moved both by sliding and by turn- 
ing over so that they were alternately convex upward and 
downward. It was noticed that the single valves which became 
buried always did so in a convex-upward attitude, this posi- 
tion seemingly being most stable in the flume. According to 
the evidence of these tests, however, there is no reason to be- 
lieve that orientation in sedimentary rocks can be demonstrated 
by the direction of convexity of a small number of valves, 
because a sudden stilling of the flow in the flume would have 
found some valves convex upward and others downward. The 
sequence of sedimentary rocks may be strongly suggested by 
consistent orientation of the direction of convexity of a con- 
siderable group of brachiopod valves (Shrock, 1948), but the 
possibility of error increases as the number of valves diminishes. 


SUMMARY 


Some terebratuloid brachiopods are moved by currents 
in a flume which are moving no faster than natural cur- 


rents found at a depth of more than 1000 fathoms in 
the Pacific Ocean. 


The very low effective density of a shell and its enclosed 
water permit it to be moved by running water even 
though it weighs hundreds of times as much as any one 
of the motionless grains of sand over which it is rolled. 
The velocity required to initiate the motion of a shell 
depends on the sphericity (shape) of the shell as well 
as on the size and weight. The smaller the effective 
density, the more important the shape. 

Shells which are more difficult to move than the sub- 
jacent sand settle into scours which are formed in the 
sand around them. 

The movement of a shell is a function of both current 
velocity and size of the subjacent sand grains. The 
same current moves a shell on one stream bed but not 
on another, depending upon the sand size. 

The orientation of a small number of movable shells 
is essentially random in relation to current direction. 
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PLEISTOCENE CLIMATE IN 
NEW MEXICO 


LUNA B. LEOPOLD 


ABSTRACT. Recently published meridional profiles of the mean tempera- 
ture of the free atmosphere provide an opportunity to apply meteorologic 
upper air data to the problem of snowlines. The modern snowline of the 
Rocky Mountains is shown to be nearly identical to the mean level of the 
0°C. isotherm for the summer months, corroborating the well-known postu- 
late that snowlines are controlled by summer temperatures. This postulate 
is applied in speculating on the annual march of temperature in glacial 
time in New Mexico. The mean annual temperature. so derived is applied 
to the hydrologic balance of a late Pleistocene lake in the closed basin of 


Estancia, New Mexico, in an attempt to check the computation for that 
lake made by Antevs. 


E line that separates areas of snow accumulation from 
those lower areas where snow disappears in summer is the 
regional snowline (Flint, 1947). The Pleistocene snowline dis- 
cussed here is the approximation furnished by the level of cirque 
excavation. 
The altitude of the snowline of Wisconsin time was approxi- 
mately 12,000 feet in the southern Rocky Mountains and de- 


creased in height northward. Data on the elevations of Wiscon- 
sin and the present snowlines have been published by Klute 
(1928) and Louis (1926). Depending on how much snow blows 
over the mountain crest, the snowline on an individual mountain 
range is not at the same elevation on opposite sides. This was 
true in Wisconsin time as well as at present and contributes to 
the difficulties of establishing the exact position of the Wiscon- 
sin snowline. 

Certain of Klute’s (1928) snowline profiles have been re- 
plotted in figure 1. His meridional profiles for the modern and 
Pleistocene snowlines are included in figure la for Japan, cen- 
tral Europe, Mexico, and the southern Rocky Mountains. A 
profile for the southern Rocky Mountains along the 110° W. 
Long. has also been plotted from the map of Louis (1926). 

On the basis of the distribution of fossil marmots, Stearns 
(1942) suggests that the snowline of north central New Mex- 
ico might have been depressed 4500 feet during Wisconsin III 
time. He thinks that although the 2000 feet of depression sug- 
gested by Klute probably was meant to apply to the larger 
Wisconsin II ice advance, it is too small even for the smaller 
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W, advance. It will be noted from figure la that Stearns’ sug- 
gestion agrees better with Louis’ profile than does the smaller 
snowline depression proposed by Klute for the area. 

All these profiles show a comparable magnitude of snowline 
depression for different parts of the world at this latitude, a 
resemblance which has been discussed in detail by Klute (1928). 
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MERIDIONAL TEMPERATURE PROFILES IN THE FREE 
ATMOSPHERE 


The modern network of upper air meteorological stations 
(radiosonde) provides data of direct interest to the problem 
of snowlines. A profile representing average upper air condi- 
tions for summer and winter from the arctic to the equator has 
recently been published by Hess (1948). His profile extends 
along the line of 80° W. Long. through Florida and Michigan. 
From Hess’ data the profile of the freezing isotherm for sum- 
mer and winter is reproduced in figure 1b. It is necessary to 
consider the possible difference between his data at 80° W. and 
the area of interest to this paper lying near the 107° W. mer- 
idian. For this reason profiles of the freezing isotherm were 
constructed for the period 1946-48 using radiosonde data from 
Albuquerque, New Mexico, Grand Junction, Colorado, Lander, 
Wyoming, and Great Falls, Montana. The average profiles for 
summer and for the complete year are plotted in figure 1b for 
comparison with the Hess data. 

It is apparent from the figure that the summer profile for 
the 107° meridian is parallel to and slightly higher than that 
for the 80° W. line. The slight difference in elevation might be 
due in part to the different period of record used for the two 
curves. Also the high mean elevation of the land surface along 
the 107° meridian implies a high elevation heating surface and 
a consequently warm air column between the surface and the 
level of the freezing isotherm. 

The mean winter profile for the 0°C. isotherm could not 
accurately be constructed along the 107° W. meridian because 
the stations considered, which are located at moderate eleva- 
tions, have mean surface temperatures in winter near or below 
freezing. To clarify the relations between the free air tempera- 
ture profiles and the snowline curves, selected profiles from 
figures la and 1b have been replotted in figure Ic. 

Before discussing the comparison of the profiles of the freez- 
ing isotherm with that of the modern snowline, the importance 
of the precipitation factor should be emphasized. A snowline 
would be lower where precipitation is heavy than where it is 
light, and in the absence of adequate precipitation, no snowline 
would occur, regardless of temperature. 

However, as a broad generalization it can be said that at a 
given elevation, any general increase in precipitation with in- 
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creased latitude in the Rocky Mountain area is small enough 
to be nearly obscured by large local variation resulting from 
different exposure, topography, and other local features. 

It seems valid, therefore, to assume that the close corres- 
pondence of the profiles of freezing temperature in the free 
atmosphere in summer and the modern snowline, indicated by 
figure lc, results from a causal relation between the two 
factors. 

If the modern snowline were controlled by winter temperature 
alone, its elevation would not decrease with latitude northward 
from about latitude 30°, and south of that latitude it would 
bend sharply upward to be parallel to the winter freezing iso- 
therm. It seems evident, then, that the modern snowline is con- 
trolled, with respect to temperatures, by summer conditions. 
This has been asserted by Klute (1928), Antevs (1928), and 
others. Their viewpoint on this matter was based primarily on 
studies of glaciers in which it has been shown that the magni- 


tude of an ice sheet is predominantly controlled by summer 
temperature. 


ESTIMATES OF POSSIBLE TEMPERATURES IN PLEISTOCENE TIME 
IN NEW MEXICO 


It is clear that any speculation concerning meteorologic 
conditions in another geologic time is subject to grave error. 
On the other hand, certain ideas of value can be gained by 
setting up a set of hypotheses and then analyzing the relations 
which would ensue. In the present analysis these hypotheses 
will, in part, be stated quantitatively. Assignment of specific 
values to factors such as temperature should be construed as a 
quantitative statement of an assumption and not as an asser- 
tion of reality. 

The first assumption made here is that the Pleistocene snow- 
line in New Mexico was controlled by summer temperature. 
For the moment variations in precipitation will be omitted 
from consideration. It is proposed to investigate the possible 
variations in snowline by assuming variations in temperatures 
alone. 

The mean lapse rate of temperature (rate of decrease with 
elevation) measured by soundings varies only slightly from 
year to year at a particular point and only slightly between 
points in a given geographic region. Conrad (1942) has shown 
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that lapse rates are remarkably similar over large portions 
of the globe. For the approximation desired here, it will be 
assumed, therefore, that the lower troposphere above New 
Mexico had the same lapse rate in Pleistocene time as at pres- 
ent. This assumption would not apply near the edge of a large 
glacier but appears to be reasonable for an area some distance 
from the edge of the major ice sheets, as was New Mexico. 

At levels near the freezing isotherm, the Albuquerque radio- 
sonde data in July 1947 and 1948 showed a decrease of tem- 
perature with height amounting to about 8 Centigrade degrees 
per 1000 meters. This rate is somewhat higher than the aver- 
age for winter at the same station, and being based on a short 
period of record, may be too high to represent an average sum- 
mer condition. To be in accord with the data for a large num- 
ber of areas studied by Conrad, a value of 6 C.° (10.8F.) per 
1000 meters will be used in the analysis. 

The second assumption is, then, that if the snowline were 
lowered 1000 meters, the July mean temperature at a given 
place would be decreased 6 C.° (10.8F.). 

It is probable that glacial times were more different from the 
present in summer than in winter. A blanket of snow over Can- 
ada and northern United States provides cooling by radiation 
just as intense as if there were a great thickness of ice under 
the snow. Though the thickness of the ice has certain meteor- 
ologic effects, particularly owing to the considerable mass of 
air which it displaces, the main difference in winter probably 
lies in the permanence of the snow cover during glacial times 
compared with the alternation of snow cover and lack of snow 
which characterizes the modern winter in the latitude band from 
40° to 50° N. Greater duration of cooling would make winter 
considerably colder over a large glacial mass and at its margin 
than at some distance south of the ice sheet. As far south as 
New Mexico, it is probable that the coldest month was little 
if any colder than the present. Thus it is postulated for pur- 
poses of discussion that the lowering of summer temperature 
was much greater than the lowering of winter temperature. 

To put this postulate in quantitative terms, it will be as- 
sumed initially that in glacial times, Santa Fe, elevation 7000 
feet, had the same average January temperature as at present. 
Later, various amounts of lowering of winter temperature will 
be discussed. Consideration of summer lapse rates leads to the 
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hypothesis that a lowering of the snowline would be accom- 
panied by a decrease of July temperature by 10.8 F.° per 1000 
meters depression. It will be assumed, further, that the mean 
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temperature of each other month would be decreased by a 
constant percentage of the difference between the temperature 
of that month and January temperature. In other words, there 
would be a graduated reduction of the temperature of each 
month, with the maximum reduction in July and no reduction in 
January. 

Such a reduction provides a new annual march of tempera- 
ture and, therefore, a new annual mean temperature. Assumin 
a snowline depression of 1500 meters or a reduction of July 
average temperature of 16.2 F.°, the new temperature curves 
for Santa Fe, New Mexico, and Fraser, Colorado, have been 
computed and are presented in figure 2a. 

Santa Fe lies near the foothills of the Sangre de Cristo 
Mountains and at the head of a broad sloping pediment surface. 
It can be considered a slope station. Fraser is in a valley nearly 
surrounded by mountains. Neither of these locations is clima- 
tologically similar to exposed mountain flank locations where 
glaciers would form, but temperature records from the latter 
locations are not available in the area. 

At Fraser, the elevation of which is 8670 feet, the present 
mean annual temperature is 32.0°F. The annual mean com- 
puted under the assumptions stated above is 24°. A uniform 
decrease of 8° in each month would have resulted in four months 
with mean temperatures well above freezing as can be seen by 
inspection of figure 2a. In modern times the same four months 
have above freezing mean temperatures. The graduated reduc- 
tion here postulated provides three months with mean tempera- 
tures less than 37° but above freezing, and the other nine 
months with subfreezing averages. The efficacy of such a 
reduction in bringing on glacial conditions seems clear. 

At the latitude of Fraser, the Pleistocene snowline lay about 
at the elevation of this station at the maximum depression. Had 
the temperature reduction amounted to 8 F.° in each month, it 
is nearly certain that a glacial condition due to temperature 
would not have been attained, while the graduated reduction 
postulated probably meets the requirement. 

Under the set of assumptions used here there would be a 
definite relation between the lowering of the snowline and the 
corresponding reduction of summer and annual temperatures. 
From the graph of figure 2b it can be seen that if the snowline 
were lowered 925 meters as postulated by Antevs (1935) there 
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would be a reduction of July temperature by about 10°F. and 
of mean annual temperature by 5 F.°. Antevs estimated the 
lowering of mean annual temperature to be 8 F.° on the basis 
of the latitudinal gradient of elevation of the modern snowline. 

It can be concluded, then, that the application of the method 
used here gives results of the same order of magnitude as those 
of Antevs and extends the reasoning to define the annual distri- 
bution of temperature. 

It is clear, however, that the determination of the tempera- 
ture of a glacial period depends directly on the value used for 
the mean depression of the snowline. Application of the present 
method to the various estimates of the snowline depression 
mentioned previously gives the following results: 


1 


Use of Various Estimates of Pleistocene 
Snowline Depression to Estimate the Reduction 
of Temperature 


Author Average Estimated Glacial 
Depression of Reduction of Mean Stage 
Snowline Temperature (F.°) 

(meters) (feet) Annual Summer 
Klute 650 2000 8.5 7.0 Ww, ? 
Louis 1800 5500 9.8 19.6 ? 
Antevs 925 2800 5.0 10.0 ? 
Stearns 1500 4500 8.2 16.4 Ws 


LAKE ESTANCIA 


In north-central New Mexico, bordered by the Sandia-Man- 
zano Mountains on the west and Mesa Jumanos on the south, 
is a broad closed basin at an elevation of about 6000 feet. This 
flat grassy valley is at present partly covered with small farms 
where the raising of frijoles (pinto beans) and small herds of 
cattle supports a scattered population. It was early recog- 
nized by Meinzer (1911) that this basin was once the bed of a 
lake, which at its maximum extent was 150 feet deep and had 
a surface area of 450 square miles. Bryan and McCann (un- 
published) have discussed in detail the origin of the basin and 
present evidence which indicates that there were actually two 
lake stages, of which the second was smaller and contained 
within the one outlined by Meinzer. 
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Small glaciers existed in Wisconsin time in the highest parts 
of the Sangre de Cristo range 100 miles north of the Estancia 
Valley. The peaks of these mountains exceed 12,000 feet in ele- 
vation, but because of their relatively low latitude, the largest 
glacier located above the present site of Santa Fe, was only 7 
miles long. Melt waters from these glaciers fed either the Rio 
Grande to the west or the Pecos River system to the east, but 
did not contribute water tothe closed basin of Lake Estancia. 

The time relations between Lake Estancia and the glacial 
chronology have been studied by Antevs (1935, 1949) in an 
attempt to date the artifacts and bones found on the Llano 
Estacado. It is one purpose of the present discussion to review 
the Antevs analysis. 

Hydrologic Balance of Lake Estancia.—Having made esti- 
mates of the relation between changes in temperature and the 
corresponding amounts of depression of the snowline in Pleisto- 
cene time, it is possible to examine the relation of these tempera- 
ture changes to the water supply of Lake Estancia. The plan 
will be as follows. Given the depression of the snowline as de- 
termined from geologic evidence, the corresponding tempera- 
ture changes will be used to estimate the change in annual 
evaporation. The resulting value of evaporation will be applied 
to Lake Estancia at its maximum stage to determine whether 
the lake could have existed under conditions having the same 
precipitation as that now experienced in the locality. If the 
lake could not be maintained under the assumed conditions of 
reduced evaporation, then it must be postulated that a greater 
rainfall existed at that time than at present, and some estimate 
of the magnitude of this increase might be made. 

Values of evaporation from open lakes in the western United 
States are meagre. Although there are a reasonable number of 
evaporation data for Class A pans, a coefficient must be chosen 
to convert even roughly the pan data into evaporation from the 
free water surfaces of lakes. 

Though there is no direct relation between air temperature 
and evaporation, a plot of mean monthly temperature against 
the mean evaporation im respective months is a graph widely 
used in hydrologic practice. The graph shows a hysteresis loop 
proceeding through the seasons. Figure 3a shows these loops 
for a number of lakes in the Great Basin determined by meas- 
urement of inflow and outflow (Harding, 1945). As might 
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be expected, evaporation generally increases with mean air 
temperature. 

It is necessary to estimate such a hysteresis loop for central 
New Mexico. No evaporation record is available for Estancia. 
The nearest record is a Class A evaporation pan at Santa Fe 
which has been plotted as the full line in figure 3b. The rain- 
fall and temperature conditions at Santa Fe are nearly iden- 
tical to those of the Estancia Valley. Though the topography 
is somewhat different, Estancia being more exposed and for 
this reason possibly more windy, the Santa Fe curve will be 
applied after adjustment. Using a coefficient of .70 to reduce 
the pan observations to conditions of a free water surface, the 
evaporation measured for each month at Santa Fe was reduced 
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to provide the hysteresis loop plotted as the dashed line in 
figure 3b. 

Using the method outlined, the annual evaporation from a 
lake surface at Santa Fe is estimated to be 44 inches. Inter- 
polating from the map of Furness (1947) a value of 48 inches 
would be obtained. His map of reservoir evaporation was de- 
rived from pan records similarly adjusted by a coefficient of 
.70. The map also indicates that Estancia at present would 
experience an annual evaporation about 3 inches greater than 
Santa Fe. Since this is much smaller than other possible errors 
in the present analysis, the Santa Fe record will be assumed 
applicable to Estancia. 

A quantitative method has been specified for computing a 
new mean temperature for each month of the year for any given 
amount of snowline depression. Curves of annual march of 
temperature computed in this manner were presented in figure 
2. In the next step, a given snowline depression was assumed 
and a new mean temperature was computed for each month. 
These values were applied to the dashed curve of figure 3b to 
obtain a new set of monthly evaporation values. By adding the 
twelve monthly values, a new annual evaporation was obtained. 
This process was performed for a number of values of snowline 
depression and a graph relating annual evaporation to the 
depression of the snowline was constructed. 

Figure 4 is a nomogram or composite graph of which quad- 
rant A is reproduced from figure 2b. It shows the relation of 
the lowering of the snowline (abscissa) to the lowering of mean 
annual temperature (ordinate). The same ordinate can now 
be applied to quadrant B where the evaporation values just dis- 
cussed are plotted. 

Present conditions are indicated by the origin of the graph 
in quadrant A; that is, zero snowline depression and zero 
temperature reduction. In quadrant B, y equals 0 when 2, the 
evaporation, is 44 inches. This value is the present annual 
evaporation of an open body of water at Santa Fe and was 
computed by adding the values around the hysteresis loop of 
the dashed line of figure 3b. 

For a value of 1500 meters of snowline depression, one can 
follow the arrows up to the sloping line of quadrant A, and 
left to the value of 8° representing the corresponding lowering 
of mean annual temperature. Proceeding horizontally from 
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y = 8° to the evaporation line and thence downward, one can 
read the value of 29 inches for the annual evaporation corres- 
ponding to a snowline depression of 1500 meters. 

Let us return for a moment to the assumption that in Pleis- 
tocene time the mean January temperature remained equal 
to that of the present time while the summer temperature was 
considerably lower. Though winter might have been somewhat 
colder than the present, for the reasons already stated it is 
probable that the reduction in winter was small compared with 
that of summer. To test the effect of a reduced winter tempera- 
ture, assume that the mean January temperature was 5F.° 
lower than at present and that the July temperature was re- 
duced 16.2F.° in accordance with a snowline depression of 1500 
meters. Under this set of assumptions the mean annual tem- 
perature would be lowered 11.8F.°. It can be seen that the new 
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assumption of a specific lowering of January temperature 
would increase the slope of the curve shown in quadrant A of 
figure 4. 

Using the new annual march curve of temperature a new 
value of annual evaporation can be obtained. The assumptions 
and results are summarized below: 

Assume: 

1) Depression of snowline of 1500 meters 

Therefore: July temperature lowered 16.2F.° 

2) January temperature lowered 5F.° and intervening 

months proportionately 

Compute from above assumptions: 

1) Annual mean temperature reduction of 11.8° F. 

2) Annual! evaporation from lake, 26 inches 

As would be expected, the additional assumption of a lower 
winter temperature would reduce the mean annual temperature. 
The consequent reduction of evaporation would be relatively 
small because the evaporation occurs primarily in summer. In 
summary, an assumed reduction of winter temperature provides 
a steeper slope to the line in quadrant A but the relations in 
quadrant B remain essentially unchanged, and for purposes 
at hand, the original assumptions provide a satisfactory 
approximation. 

As Lake Estancia lay in an enclosed basin, the rainfall on 
the lake plus runoff into the lake was balanced by the evapora- 
tion from the lake surface. 

Let e equal annual evaporation in inches 
r equal annual runoff in inches 
p equal annual precipitation 

The contributing drainage area was 1550 square miles, and 
the area of the lake at maximum stage was 450 square miles. 
Following the reasoning used by Antevs (1935), the total vol- 
ume of water can be equated: 


p (450) + r (1550) = e (450) 


P\ 1550) 


This is a linear equation the solution of which may be repre- 
sented as a family of straight lines on a graph. It is shown in 
quadrant C of figure 4. Given values of evaporation (abscissa) 
and precipitation (ordinate), the runoff r necessary to main- 
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tain the lake can be read off the graph by interpolation between 
the sloping lines. 

It is obvious that for any given value of annual precipitation 
only a certain range of annual runoff values can occur in na- 
ture. A map of annual runoff for the United States has been 
published by Langbein (1949). It shows values of .25 inches 
near Estancia, .50 inches near Santa Fe, and an increase of 
values with elevation in central New Mexico. The relation of 
rainfall to runoff existing under present vegetative types and 
the geologic formations occurring in New Mexico determines 
the possible range of values, and this range is indicated by the 
stippled area in quadrant C of figure 4. 

It is now possible to determine the rainfall necessary to 
maintain Lake Estancia under various assumed conditions of 
evaporation which in turn are determined by temperature. Re- 
ferring back to the example cited above, a depression of the 
snowline of 1500 meters, one can follow the arrows along the 
dashed line into quadrant C and observe that according to the 
assumptions made, a precipitation of 21 inches would be neces- 
sary to maintain the lake at its high stage. The reasoning so 
far has considered no change of precipitation but only changes 
of temperature. It can be reasoned, then, that the snowline de- 
pression of 1500 meters indicated by geologic data and caused 
by temperature alone would not be associated with a low 
enough evaporation to maintain Lake Estancia under the pres- 
ent precipitation of 14 inches. If the precipitation were 21 
inches at that time, the increased precipitation would also ap- 
ply in the mountains and would in itself lower the snowline. 

Therefore, the case of no temperature reduction will be con- 
sidered. If the temperature were the same as that of the present, 
the lake evaporation would be 44 inches (y = 0, « = 44, in 
quadrant B). To maintain Lake Estancia under that condition, 
the precipitation would have to be about 30 inches. Such a 
rainfall would have brought the upper part of the pine forest, 
bordering the spruce zone, down to the margins of Lake Estan- 
cia. Such an occurrence would surely be indicated by well-de- 
veloped soil profiles, which are not to be found anywhere in the 
Estancia area. Actually, some streams tributary to Lake Es- 
tancia have not discharged enough water since the lake receded 
to breach the lake shore deposits damming them (Bryan and 
McCann, manuscript). It is impossible to believe, then, that 
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increased precipitation without, lowered temperature could ac- 
count for the lake. 

Great variations in rainfall would change the vegetative as- 
sociations sufficiently to affect the runoff coefficient. Also the 
change in the period of snow cover might have an important 
effect on runoff. However, as long as the postulated changes 
require only moderate changes in precipitation and do not re- 
quire major changes in the vegetation in the Estancia Valley, 
the use of modern rainfall-runoff relation in the present analy- 
sis seems justified. 

The analysis leads to the conclusion, therefore, that Lake 
Estancia must have existed at a time when temperature was 
lower and precipitation higher than the present. How much 
did each of these factors contribute to the lowering of the snow- 
line? A postulation of the relation of temperature, explored in 
some detail in the present study, has been made possible owing 
to two basic facts: first, the modern snowline seems controlled 
to a considerable degree by summer temperature, and second, 
tht average lapse rate of temperature is so nearly uniform 
around the earth in the middle latitudes that an assumption of 
similar lapse rate in Pleistocene time seems very reasonable. 

The relation of snowline to precipitation is much more dif- 
ficult to postulate because it has not been investigated in detail 
in the Rocky Mountains even for modern conditions. Although 
a reasonable guess might be made, it would be founded on an 
even less logical set of assumptions than the temperature rela- 
tion already discussed. 

Consider the graphs of figure 4 for an assumed depression of 
the snowline of 1000 meters. Following the dashed line around 
the graph from that value, it is seen that a precipitation of 
24 inches is necessary for maintaining the lake. That the dif- 
ference between 24 inches and present 14 inches would account 
for the extra 500 meters of snowline depression necessary to 
make the 1500 meters inferred from geologic evidence, is a 
speculation which cannot be supported or tested. 

It seems desirable to mention one other feature of the present 
analysis. Antevs (1935) used very similar reasoning to that 
presented here and concluded that Lake Estancia could not 
have existed at the same time that glaciers occurred in the 
mountains. He postulated (Antevs, 1949) that the lake post- 
dated the maximum glaciation and existed 10,000 to 12,000 
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years ago at a time when both temperature and precipitation 
were higher than the present. He chose values of 20 inches for 
precipitation and 35 inches for evaporation as likely during the 
existence of the lake. These values require a runoff of about 4.5 
inches from a precipitation of 20 inches, a high runoff co- 
efficient for the western country. Assuming, however, that this 
might still be possible, his postulation of a temperature higher 
than that of the present would necessitate an increase in evap- 
oration. The evaporation data from Santa Fe used in the pres- 
ent study indicate that 35 inches of evaporation would require 
a lower, not a higher temperature. Antevs used the evaporation 
of Lake Tahoe, California, as a criterion. The question is 
whether the Lake Tahoe evaporation provides a better measure 
of evaporation in New Mexico than an adjusted Santa Fe 
record. 

The present analysis supports the conclusion of Antevs that 
Lake Estancia must have existed in a period more humid than 
the present. Antevs interpreted the necessity for increased pre- 
cipitation as an indication of a post-glacial date, while the 
present study leads to the conclusion that a decreased tempera- 
ture should also be necessary, indicating contemporaneity of 
the lake with glaciers in the mountains. These opposing view- 
points point to the hazard of using a date assigned to the lake 
as a measure of the age of associated artifacts. 

The writer gratefully acknowledges the help of Prof. Kirk 
Bryan, Prof. Chas. F. Brooks, and Dr. V. Conrad, of Harvard 
University, and Mr. G. E. Harbeck of the U. S. Geological 


Survey, who read the manuscript and made many valuable 
suggestions. 
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COMMUNICATIONS 


LYELL, GREEN, HILGARD, AND THE 
MISSISSIPPI VALLEY LOESS 


Sir Charles Lyell published the tenth edition of his Principles of 
Geology at London in 1867. The first edition had been published 
in 1880-82. On Sept. 25, 1866, Sir Charles had written to Col. 
James Green of Dayton, Ohio, the holograph letter quoted below. 
This letter, with the original envelope and its postage stamp, are 
pasted in the front of the first of the two volumes of the tenth 
edition of Lyell’s Principles which now belong to the present writer. 

“Dear Sir,” Sir Charles wrote to Colonel Green from 53 Harley 
Street, London, “I have just received your letter informing me of 
the very interesting discovery of the remains of the buffalo fish in 
the loess. Although the chapter in which I shall be able to allude 
to it, is in the hands of the printer it will not be too late to get in a 
reference to it. I am glad you are pursuing your researches in that 
direction. It will throw great light on the former geological history 
of the Valley of the Mississippi—As my proofs are coming in fast 
& I am keeping up an active correspondence, you will oblige me 
much by letting Professor Hilgard of Oxford Mississippi know 
that I have received his valuable letter a few days ago. I am sorr7 
that he did not know that my book has taken a longer time to get 
through the press than I had anticipated as I should like very much 
to have iearnt from Mr. Conrad whether the only two species from 
the Artesian well of New Orleans which he (M*‘ Hilgard) had 
been unable to identify were known as living shells or not. They 
were I think a Cardium & Tellina. I fully expect that they will 
turn out to be recent species because the Gnathodon cuneatus & 
several other living shells are found at a much greater depth Even if 
they should not be known shells they might be living in the Gulf 
for we have found Murese Vaginatus which I have figured in Chap 
XIII of the last edition of my Elements in the Mediterranean. It 
is now excessively rare but was evidently abundant towards the 
close of the Newer Pliocene Period. 

“But I must not indulge in this subject at present but hope to 
write to M* Hilgard at some later day & in the mean time I beg to 
thank him for his prompitude in sending me the latest information 
which he had obtained. 

“IT am, dear Sir 

“very truly yours 
“Cha Lyell” 

Colonel Green added an asterisk after the word “obtained” in 
Sir Charles Lyell’s final paragraph, and then wnote on the opposite 
page, below Lyell’s signature, the following note. 
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“P.S. It was Prof Hilgard who urged me to forward immediately 
to Sir Charles Lyell, a Description of the fossil remains that I 
discovered in the Loess of the Mississippi Near Vicksburg so that 
It might be published in the New addition of his ‘Principles of 
Geology’.” 

“James Green” 


Upon a number of pages of my copy of the 1867 edition of Lyell’s 
Principles of Geology there are interesting pencilled annotations by 
Colonel Green, for example on page 464 of Vol. I where the buffalo 
fish in the loess near Vicksburg are described by Lyell. It may 
not be without interest to observe that the 1867 copy of Lyell’s 
Principles which contains the holograph letter quoted above is also 
annotated with the words “Col. James Green with the compliments 
of Wm. M. Smith purchased in London. Aug. 1868.” One wonders 
whether still other evidences of early correspondence between 
American geologists and European geologists are extant. For Lyell 
letters to Rev. Charles Kingsley and to Rev. Charles Pye Smith, 
both Englishmen, see Henry P. Zuidema, Jour. Geology, vol. 55, 
no. 5, Sept., 1947, pp. 439-445. 

LAWRENCE MARTIN 


Wasnineton, D. C. 
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ENROLLMENT OF GEOLOGISTS AND GEOPHYSICISTS 
IN THE NATIONAL SCIENTIFIC REGISTER 


Planning for the mobilization and effective use of the scientists 
and engineers of the United States for purposes of national defense 
is now under way. Appointment of a committee composed of leaders 
in science and industry by the National Security Resources Board 
is a recent step. The representative of the geological professions 
on this committee is E. De Golyer; the chairman is Charles A. 
Thomas, vice-president of Monsanto Chemical Co. and former 
president of the American Chemical Society. 

A basic element in the entire program for maximum use of the 
nation’s scientific manpower will be an adequate scientific register. 
This will tell us the magnitude and nature of the national pool 
of technological competence. 

In order to secure this information, the N.S.R.B. was instrumental 
last summer in causing the establishment of a National Scientific 
Register Office in the U. S. Office of Education. This office has let 
contracts with the American Chemical Society, the National Research 
Council and the American Institute of Physics for enrollment of 
all physical and biological scientists and mathematicians. Enroll- 
ment will be accomplished through the sccieties and the register 
will also be kept up to date through the societies. The American 
Geological Institute has been asked to ensure the enrollment of all 
geologists and geophysicists. 

The project must start essentially from scratch because the 
records of the World War II National roster are completely out of 
date. Also, with the tremendous increase in the number of working 
scientists since 1945, its records are very incomplete. 

Some geologists and geophysicists who are included in American 
Men of Science, filled out a comprehensive registration form in 
1948. These men will not be asked to answer all the questions on the 
new simpler form, but will be asked to fill out basic questions in 
the interest of uniformity. 

The following methodology is being established for this project 
to ensure effective coverage. It is hoped that use of this method will 
spread the work, prevent wasteful and annoying duplication, and 
get the job done more promptly and effectively. 

The register will include all geologists and geophysicists with a 
bachelor’s degree (or equivalent) and all graduate students. 

Circularization will be accomplished by district committees. 
These will center in active local societies, state geological surveys 
and major departments. 


The register will be used for the following purposes— 
(a) by Scientific Advisory Committees in the Selective Service 
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System. It is expected that these committees will review appeals for 
deferment which are made to The President and that the informa- 
tion on file with the Register will be used in the consideration of 
such cases. 

(b) by the National Security Resources Board. This agency is 
charged with responsibility for the planning of general mobilization 
both military and civilian, and will use the information in determin- 
ing the requirements for scientific manpower in various activities. 

(c) by the National Research Council and by the American 
Geological Institute in determining the distribution and availability 
of earth scientists among the various specialized occupations. 

(d) as a means whereby Government agencies may identify per- 
sonnel with specialized training or experience for unusual work 
requirements. 

(e) in the preparation of studies and recommendations to Govern- 
ment agencies, such as the National Security Resources Board, 
concerning the needs for geological scientists in civilian occupations. 

In addition, it will furnish information concerning the composi- 
tion and structure of the geological professions which will be 
invaluable. 

The full cooperation of every member of the geological professions 
is earnestly requested so that our part in this national effort may be 
performed quickly and effectively. 


‘ 
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Elementary Pile Theory; by H. Soopax, and E, C. Campse.i. 
Pp. ix, 78. New York, 1950 (John Wiley & Sons, Inc., $2.50). 
—This booklet is noteworthy in two important respects. First, it 
opens to physicists at large a field still somewhat shrouded in 
mystery, but naturally interesting to students who are alive to 
modern research, And it may be said, to the credit of the authors, 
that their treatment will doubtless attract readers by its clarity, 
directness and simplicity. Secondly, the little volume revives a 
part of classical physics that is frequently neglected in the training 
of physicists during these times of formalized pursuits. For it 
presents in palatable fashion the elements of the kinetic theory, 
bringing out the visual beauty of its concepts without the dry appeal 
to transport equations, distributions in phase space and the like. 
The contents of the book, if injected discriminately into a standard 
course in the kinetic theory, would enliven it greatly. 

Using fairly simple mathematics (i.e. nothing beyond calculus) 
the authors treat the essential details of the functioning of chain- 
reacting piles. Illustrations are used freely in explanation of the 
physical and mathematical ideas encountered. The style is precise 
and non-discursive, but the economy of statement is not carried to 
the extreme which often discourages readers of technical treatises. 
Still, if the book has any fault at all, it might lie in its confinement 
to issues directly related to the pile, in its failure at times to point 
to other situations in which similar reasoning is applicable. But the 
brevity of the text makes this restriction unavoidable. 

HENRY MARGENAU 


Stratigraphy and Paleontology of the Brownsport Formation 
(Silurian) of Western Tennessee; by Tuomas W. Amspen. Pea- 
body Museum of Natural History, Yale University, Bulletin 5. 
Pp. 188; 80 figs., 84 plates. New Haven, 1949 (Yale University 
Press, $5.00).—The Brownsport formation has long been a center 
of paleontologic interest as a collecting ground for ‘fine Silurian 
fossils. Since the pioneer work of Troost, many geologists have 
added to our knowledge of Brownsport rocks and fossils, notably 
Roemer, Safford, Foerste, Pate and Bassler, and Springer. Amsden’s 
report presents the results of a detailed stratigraphic and paleonto- 
logic study. An outstanding contribution is the careful description 
and illustration of 97 species of brachiopods, corals, cystoids, and 
blastoids, of which 87 species and one genus are new. Except for 
Springer’s work on the crinoids, the Brownsport fauna has never 
been systematically studied, so that with the publication of Amsden’s 
report a description of this rich and varied fauna is available for 
the first time. With this as a basis, Amsden documents the correla- 
tion of the Brownsport formation with the Bainbridge limestone 
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of Missouri, the Louisville limestone of Kentucky, the Henryhouse 
shale of Oklahoma, and part of the Niagaran dolomites in the 
Chicago area. 

In its stratigraphic conclusions this report differs from the 
interpretation first made by Pate and Bassler in 1908, and later 
endorsed by other workers, that the Brownsport is to be considered 
a group comprising in ascending order the “Beech River formation,” 
“Bob limestone,” and “Lobelville formation.” A summary state- 
ment of the Brownsport published as part of a larger work by C. 
W. Wilson, Jr. in 1949 follows this interpretation. Amsden believes 
that the basis of this subdivison is invalid. Whereas the “Bob” has 
been considered as a sheet of massive limestone separating the thin- 
bedded, argillaceous limestones and shales of the “Beech River” 
and “Lobelville” formations, Amsden’s work indicates that the 
“Bob” comprises rather local and lenticular facies equivalents of 
“Beech River” and “Lobelville” lithology. This interpretation of 
the stratigraphy, based on detailed examination and measurement 
of 53 sections, together with a biostratigraphic analysis indicating 
that Bassler’s formations are not faunally distinct, prompts Amsden 
to abandon the tripartite division. JOHN IMBRIE 


The Geology, Paragenesis, and Reserves of the Ores of Lead and 
Zine; Report of the 18th Session of the International Geological 
Congress, part 7, Great Britain, 1948; K. C. Dunnam, Editor. Pp. 
400, 99 figs. London, 1950.—The International Geological Con- 
gress has made another valuable contribution to economic geology 
in compiling a wealth of information on the geology, paragenesis 
and reserves of the ores of lead and zinc. 

The urgent need for such a report is well pointed out by the 
opening remarks of the chairman, Professor W. R. Jones: “Half of 
the world’s supply of lead and zinc is produced in North America 
.. . It is over a quarter of a century since the last discovery of a 
lead and zinc deposit of major importance . . . Only a few of 
the major fields contain proved reserves sufficient at their present 
rate of production to last twenty years or so . . . Production dur- 
ing the last twenty-five years has been greater than during all 
previous mining history . . . Here, then, is a challenge to the 
geologists and mining engineers of the world to find more deposits 
and to increase production.” The object of the symposium is to 
“help the geologists meet this challenge by increasing our knowledge 
concerning the geological conditions most suitable for the occurrence 
of the different types of these ores.” 

Lead and zine deposits “have been formed in the earth’s crust at 
levels varying from some of the deepest which can be observed up 
to the top of the permanent groundwaters.” A fivefold classifica- 
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tion is given—hypothermal in regional-metamorphic environments ; 
pyrometasomatic-mesothermal; deep mesothermal veins and re- 
placements; leptothermal veins and replacements; telethermal re- 
placements and disseminations. Hypothermal deposits contain the 
largest reserves. 

A very comprehensive table has been prepared of the principal 
lead and zinc mines of the world, their country rock, associated 
igneous rocks, mineral paragenesis, nature of ore shoots, reserves, 
assays of ore and current production. Although the production and 
reserve figures are subject to the limitations imposed by a de- 
ficiency of information, the table is a valuable summary of the 
world’s literature. 

The report is presented in place of the usual volume on world 
resources, as prepared for other sessions. However, it follows the 
pattern of the earlier reports in presenting well-written descriptions 
of individual lead and zine deposits. Each description is written by 
authors well qualified for the task, and includes an abstract in 
English and a bibliography. 

Most countries which hold important positions as lead and zinc 
producers have contributed papers. Russia is a notable exception. 
Forty-four authors have contributed 83 papers which may be 
divided as follows: 2 from Canada, 7 from the United States, 2 
from Mexico, 1 from Peru, 2 from Australia, 1 each from Tangan- 
yika, Nigeria, Morocco and Tunisia, 2 from Portugal, 4 from Italy, 
1 each from Jugoslavia, Greece and Poland, 2 from Germany, 1 
from Great Britain, 2 from Sweden and 1 from China. Papers on 
newly described districts are presented, including the Taxco Mining 
District, Mexico, and deposits in French Morocco, Tunisia and 
Tanganyika. New information on the Mississippi Valley, Franklin 
Furnace, Bisbee and others bring our storehouse of knowledge up to 
date. An excellent description of the Cerro de Pasco properties is 
presented. The Kupferschiefer is described and a supergene origin 
reafirmed. Several authors suggest that supergene leaching of lead 
has taken place. A tendency to retain theories of supergene deposi- 
tion of sulfides is shown by the statement of one author “ .. . no 
igneous rocks of whatever kind can be found. The ore consists 
chiefly of galena in which silver content is usually high. Gangue 
minerals are calcite and quartz. This type of deposit represents, 
probably, the product of groundwater concentration of non-igne- 
ous origin.” Theories on sedimentary lead and zinc deposits seem 
to be restricted to Europe. 

The publication represents an international exchange of facts and 
theories which will work to the benefit of all geology, and espe- 
cially those geologists directly concerned with lead and zinc. 

CHARLES H. SMITH 
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PUBLICATIONS RECENTLY RECEIVED 


Elektronenoptik I. Grundstige der Theoretischen Elektronenoptik; by Alex- 
ander A. Rusterholz. Basel, 1950 (Verlag Birkiiuser, paper, 25 Swiss 
francs; cloth bound, 29 Swiss francs). 

Illinois State Geological Survey: Report of Investigations No. 146, Appli- 
cation of Mineralogy to Soil Mechanics; by Ralph E. Grim; Circular No. 
149, Where to Find Information on Mineral Raw Materials; by M. M. 
Leighton; Circular No. 156, Acid Etching in the Study of Limestones 
and Dolomites; by J. E. Lamar; Urbana, 1950. 

Grassland Historical Studies: Natural Resources Utilization in a Back- 
ground of Science and Technology. Vol. 1, Geology and Geography; by 
James C. Malin. Lawrence, Kansas, 1950 (printed by the author, 1041 
University Drive, paper cover, $2.50). 

Applied Sedimentation; edited by Parker D. Trask. New York, 1950 (John 
Wiley & Sons, Inc., $5.00). 

Thermodynamics. An Advanced Treatment for Chemists and Physicists; 
by E. A. Guggenheim. Monographs on Theoretical and Applied Physics, 
Vol. Il. New York and London, 1950 (Interscience Publishers, Inc., 
$6.50). 

Experimental Designs; by W. G. Cochran and G. M. Cox. New York, 1950 
(John Wiley & Sons, Inc., $5.75). 

Introduction to Theoretical Igneous Petrology; by Ernest E. Wahlstrom. 
New York, 1950 (John Wiley & Sons, Inc., $6.00). 

Theorie und Lésungsmethoden des Mehrteilchenproblems der Wellenme- 
chanik; by P. Gombfs. Basel, 1950 (Verlag Birkhauser, paper, 24.50 
Swiss francs; cloth bound, 29.50 Swiss francs). 

The Nature of Natural History; by Marston Bates. New York, 1950 
(Charles Scribner’s Sons, $3.50). 

Early Man in the New World; by Kenneth Macgowan. New York, 1950 
(The Macmillan Company, $5.00). 

Die Lebensweise der Dinosaurier; by M. Wilfarth. Stuttgart, 1949 (E. 
Schweizerbart’sche Verlagsbuchhandlung, paper bound, D. M. 12), 

Geochemistry; by Kalervo Rankama and Th. G. Sahama. Chicago, 1950 
(The University of Chicago Press, $15.00). 

Measuring Our Universe; by Oliver J. Lee. New York, 1950 (The Ronald 
Press Company, $3.00). 

The Nature of Physical Realty; by Henry Margenau. New York, 1950 
(McGraw-Hill Book Company, $6.50). 

This Earth of Ours—Past and Present; by C. Wroe Wolfe. Revere, 
Massachusetts, 1950 (The Earth Science Publishing Company, $5.00). 
Origin of Kansas Great Plains Depressions; by John C. Frye. Kansas 

State Geological Survey Bulletin 85, Part 1. Lawrence, 1950. 

Introducing the Insect; by F. A. Urquhart. New York, 1950 (Henry Holt 
and Company, Inc., $5.00). 

The Chemical Elements and Their Compounds, 2 vols.; by N. V. Sidg- 
wick. Oxford and New York, 1950 (Oxford University Press, $14.00). 

First Principles of Atomic Physies; by Richard F. Humphreys and Robert 
Beringer. New York, 1950 (Harper & Brothers, $4.50). 

Textbook of Biochemistry, 5th ed.; by Benjamin Harrow, Philadelphia and 
London, 1950 (W. B. Saunders Company, $6.00). 

Laboratory Manual of Biochemistry, $d ed.; by Benjamin Harrow, G. C. 
H. Stone, Harry Wagreich, Ernest Borek, and Abraham Mazur. Phila- 
delphia and London, 1950 (W. B. Saunders Company, 2.25). 

Fundamentals in Chemistry for the Laboratory, revised ed.; by W. C. 
Fernelius, A. B. Garrett, and L. L. Quill. Boston, 1950 (Ginn and Com- 
pany, $2.75 paper bound). 
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